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Preface 



These lecture notes are written as reference material for the Advanced Course "Hydrodynamical 
Methods in Last Passage Percolation Models" , given at the 28*^ Coloquio Brasileiro de Matemdtica 
at IMPA, Rio de Janeiro, July 2011. 

The notes focus on recent progress made in two dimensional Poisson Last Passage Percolation 
(LPP), using what is generally known as "hydrodynamical methods" [U O [U [H [lOl [11]. These 
methods have been developed in recent years by the authors, but also by many others, such as 
Groeneboom, Seppalainen, Balazs, Aldous, Diaconis, Martin and Ferrari. The methods try to use 
probabilistic tools to analyze last passage percolation, such as Markov processes, interacting parti- 
cle processes, concentration inequalities etc. This is somewhat in contrast to the analytic methods 
used by Baik, Deift and Johansson in their ground breaking paper [2j, which were further developed 
by many other authors as well, such as for example Rains, Praehofer, Spohn and O'Connell. The 
analytic methods (which we will not discuss in any detail) are able to get strong results on limiting 
distributions and fluctuation probabilities of the length of a longest path in the characteristic di- 
rection for some specific LPP models, using exact formulas for certain distributions and analyzing 
the limiting behavior of rescaled versions of these formulas. The hydrodynamical methods are not 
able to find these limiting distributions (so far), but they were able to find the correct scaling for 
these specific models, and extend this scaling to the transversal fluctuations of the paths, giving 
information on second class particles, for example. Furthermore, it turns out that these methods 
can be (partially) extended to models for which there are no analytical methods available yet. Also, 
certain tightness properties of rescaled functionals of the LPP model can be established using the 
hydrodynamical methods, whereas limiting distributions for these functionals can be determined 
using the analytic methods. 

The main object of study is LPP using iid weights on a homogeneous Poisson process in M^. 
The strongest results will be valid for the case where all weights are equal to 1 (called the classical 
Hammersley process). Almost all results can (and have been) extended to LPP on with iid 
exponential weights, which corresponds to the Totally Asymmetric Exclusion process. However, 
for the purpose of this one week course, we will not go into this specific model. 

Chapter [1] introduces the model and shows the existence of Busemann functions, a concept 
originally developed for Riemannian geometry, but which proves to be very useful in our con- 
text. In Chapter [2] we introduce the Hammersley Interacting Fluid System, a generalization of the 
Hammersley Interacting Particle System, which turns out to be an important tool in studying the 
Poisson LPP. Chapter [3] connects the Busemann functions and the fluid process in the sense that 
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the equilibrium measures for the fluid process can be described by the Busemann functions. This 
leads to a description of the Busemann function in the classical Hammersley process and to a new 
description of the multi-class particle (or fluid) system. Chapter [4] introduces second-class particles 
in the fluid process, and shows they satisfy the usual strong law when in equilibrium, or close to 
equilibrium. It also shows how the Busemann functions can be used to describe the asymptotic 
speed of second class particles in a rarefaction fan. Finally, Chapter [5] introduces some remarkable 
identities that link the variance of the longest path to the variance of the starting configuration. 
In the classical Hammersley process, this leads to normal central limit theorems for longest paths 
in non-characteristic directions, and to cube-root fluctuations of the longest path along the char- 
acteristic direction. 
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Chapter 1 



The Poisson Last Passage Percolation 
Model 



1.1 The Graphical Construction 

Let P C be a two-dimensional Poisson random set of intensity one. On each point p G P 
we put a random positive weight ojp and we assume that {ojp : p G P} is a collection of i.i.d. 
random variables, distributed according to a distribution function F, which are also independent 
of P. When F is the Dirac distribution concentrated on 1 (each point has weight 1; we will denote 
this F by 5i), then we refer to this model as the classical Hammersley model [1]. 

For each p, q G M?, with p < q (inequality in each coordinate, p ^ q), let n(p, q) denote the 
set of all increasing (or up-right) paths, consisting of points in P, from p to q, where we exclude 
all points that share (at least) one coordinate with p. So we consider the points in the rectangle 
(p, q], where we leave out the south and the west side of the rectangle. In this probabilistic model, 
the "metric" (or last passage time) L between p < q is defined by 

L(p,q) := max { J]] Wp'} . 

roen(p,q) ^ 

Then L is superadditive, 

i(p,q) > i(p,z) +L{2,q) . 

When we consider a path vj from p to q consisting of increasing points (pi, . . . ,Pn), we will 
view zu as the lowest increasing continuous path connecting all the points, starting at p and ending 
at q, and then excluding p. This way we can talk about crossings with other paths or with lines. 
Suppose uJi and ti72 are two different paths between p and q that both have the maximal weight 
L(p,q). Now suppose that (xi,ii) < (2:2,^2) are two points on the intersection of the two paths, 
such that the paths do not intersect in the open strip (xi, X2) x M. Then one path must lie entirely 
below the other path in this strip, and the total weight of the two paths in this strip must be equal 
(otherwise we would be able to construct a "longer" path). This means that we can define the 
minimum of wi and tz72 by choosing the lowest path in each such strip. The finite geodesic between 
p and q is given by the lowest path (in the sense we just described) that attains the maximum 
in the definition of L(p,q). We will denote this geodesic by w{p,q) (this is well defined for any 
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ordered pair (p, q), even if we do not specify the order). 

1.2 The shape Theorem 

A crucial result for LPP models is the following shape theorem: set = (0,0), n = (n,n), 

F(x) = P(wp < x) and 7 = 7(F) = sup > . 

n>l n 

Theorem 1.2.1 Suppose that 

f'OO 

/ - F{x) dx<+oo. (1.2.1) 

Jo 

Then 7(F) < 00 and for all x,t > 0, as r ^ 00, 

EL{0,{rx,rt)) 



^Vxt , 
r 

and, a.s., 

L(0,(rx,rt))^^^^ ^^^.2) 
r 

Theorem 11.2.11 shows that L has a curved limiting shape, mainly due to the invariance of the 
Poisson process under volume preserving maps: if x, t, A > and p G M^, then 

L(0, (x, t)) = L (p, p + (Ax, t/A)) (1.2.3) 

(choose A = -s/t/x and p = 0). 

Proof: Equation ()1.2.3p shows that it is enough to prove the theorem for (x,t) = (1,1). When 
considering only one ray, the convergence of L is a standard consequence of Liggett 's version of the 
Subadditive Ergodic Theorem |19j . as soon as we can show that 

EL(0,(r,r)) 
lim sup < 00 . 

r— ^-oo f 

Recall that 5i denotes the Dirac distribution concentrated on 1. For each p G [0, 1] denote by 
Ep expectation for the Hammersley last passage model induced by Bernoulli weights it^p, where 
P(tt;p = 1) = p. This coincides with the classical Hammersley model, but with Poisson intensity p 
(instead of 1). Thus, it is not hard to see that 
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for some 7(5i) < Now we use an idea introduced in |20j : 

L(0,p) = max | > w^'} 

p Gro 

^en(o,p) Jo p,g^ 

/>oo 

The integrand in the last hne corresponds to the Bernoulh model with p = 1 — F(x). This means 
that 

EL(0,(r,r)) ^ Ei„^(,)L (0, (r, r)) 
hm sup < lim sup / ax 



7(5i) / Vl - F(x) . 

□ 



Theorem 1.2.2 7/ ()1.2.ip is strengthened to: there exists a > such that 

/•oo 

I exp(ax) dF(x) < +CX) , (1-2.4) 

JO 

then there exist constants oq, ai, Q^3) > such that for all r > ao 

P(|L(0, (r, r)) - 7r| > n) < ai exp f - ) (1.2.5) 

V A/riogr/ 

for n G [ 03 "v/r log^ r , a4r'^/^ log r ] . 

The proof of (jl.2.5p . which was first stated in [9], follows Kesten's approach developed for lattice 
First Passage Percolation models [18j , which is based on the Method of Bounded Increments applied 
to L (see Lemma [1.2.31 below) . This gives a bound on the fluctuation of L around its expectation. 
Then adapting a clever argument used by Howard and Newman in [T7] shows that one can replace 
EL(0, (r, r)) by the shape function. 

Lemma 1.2.3 Let {J-k}o<k<N be a filtration and let {Uk}o<k<N be a family of positive random 
variables that are Tn measurable. Let {Mfc}o<fc<Ar be a martingale with respect to {-?7c}o<fc<Af • 
Assume that the increments := Mk — Mk-i satisfy 

\^k\ < Co for some cq > (1.2.6) 

and 

E(Ai| <E(C/fc| J-fc_i) . (1.2.7) 

^It is actually known that 7((5i) = 2 [!]■ 
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Assume further that for some constants < ci, C2 < oo and xq > e^c^ we have 

fl.2i 



^ N \ 

Uk > X \ < ci exp(— C2x) when x > xq . 

\k=l J 



Then irrespective of the value of N , there exists universal constants < 03,04 < 00 that do not 
depend on iV, co,ci,C2 and xq, nor on the distribution of {Mk}o<k<N o-i^d {Uk}o<k<N , such that 

F{Mn - Mo > x) < cJl + ci + ^] exp ( -04^ ) (1.2.9) 



C2X0 ) \ 

whenever x < C2 Xq^'^ . 

Proof: See Theorem 3 in [18]. 

□ 

We decompose L as a sum of martingales increments as follows. For each integer r > 1, let 
Lr := L{0,{r,r)) and consider a partition of the two dimensional square 

[0,rf=ufL,Bi 

into N = r"^ disjoint squares of size one. Let J-q = {0, fi} and for each k = 1, . . . ,N consider the 
cr-algebra 

Jfc = a ({ (p, Wp) : p G Uti^i n P}' 

and the Doob martingale 

Mk :=E{Lr I J-fc) . 
Denote by P/ the probability law induced by 

{(P,^p) : pG5/nP}, 

and by Q,i the underlying sample space. For uj,a € YliLi 



N 



[uj,a]k := {uji, . . . ,ujk,(^k+i, ■ ■ ■ i^-rO £ JJf^r 



1=1 

Then 



Ak{uji,...,ujk) ■= Mk-Mk-i 

= J Lr[io,a]k n dFiiai 



l=k+l 

N 



Lr[uj,a]k-i'[ld¥i{ai) 
l=k 
N 

/ Lr.[uj,a]k - Lr[uj,a]k-iY[dFi{ai) 



l=k 
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(To integrate Mk over ak does not change it, and it allows us to put Mk and M^-i under the same 
integral.). 

For each 1 < A: < iV let 

Zk := ^ Wp. 

Then {^fc}i<A;<Ar is an i.i.d. collection of random variables such that 

Ee"^i = exp (Ee""' - 1) < oo . 



Lemma 1.2.4 Let 1^ denote the indicator function of the event that the geodesic Wr '■= tz7(0, (r, r)) 
has a point p £ DP . Then 

\Lr[uj,a]k - Lr[u:,a]k-i\ < ma,x{Ik[uj, a]k, Ik[i^, o-]k-i} x 

max{Zk{uJk),Zk{ak)} . 

(1.2.10) 

Proof: We note that there will be no difference between Lr[uj,cr]k and L^. [w, o"]fc_i, if no geodesic 
has a point p G -B^nP (recall that [a;, a]k and [oj, cj]fc_i only differ inside Bj.), and this corresponds 
to the first factor in the right hand side of (jl.2.10p . And, if one of them does intersect, then the 
increment can not be greater then the total weight inside the box B^. (Compare with (2.12) in 

m-) 

□ 

The next step is to construct, from {{p,Wp) : p G P}, a new process {{p,idp) : p G P}, by 
truncating the original model inside each box Bk, if Zj^ > blogr, in order to have 

Wp < b log r . 

We do this truncating by putting all the weights in box B^ equal to if the original weight is 
greater than 61ogr. Note that the truncated process in each box is still independent of all the 
other boxes. Let us denote the configuration induced by the truncated model by a) and for each 
random variable Y that depends on uj, let us write Y{uj) := Y{u}). Thus, 

\Lr[u},a]k - Lr[uj,a]k_i\ < (61og r) max {4[a;, cr]^, 4[a;, o-]fc_i } , 

and hence 

\Akiui,...,Uk)\< blogr. (1.2.11) 
The upper bound (jl.2.1ip allows us to apply Lemma ll.2.3l to get concentration inequalities for 



Mat - Mo = Lr{io) - ELr{L0) 
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Lemma 1.2.5 Let Uk := 2(6 log r)^/^. Then < 61ogr and 

E{Al I Tk-i) < E{Uk I Tk-i) . 
Proof: By Schwarz's inequality and Lemma 11.2. 41 

E(A| I J^fc-i) < / max {lk[uj,a]kJk[i^,(T]k-i} 



N 



max {Zkiojk), Zk{ak)Y'[ld¥i{<Ti) dfk {• 



l=k 



< 



{Ik[uj,a]k + Ik[i^,a]k-i) 



N 



{blogrf'[[dfi{ai)d¥k{uk) 



l=k 



K{Uk I J'k-i) ■ 



□ 



Since the number of boxes intersected the up-right path zur is 2r — 1 (with probability one), we 
conclude that 

N N 

^Uk = 2(61ogr)2^7fc < 2(61ogr)2(2r-l). (1.2.12) 

k=l k=l 

By choosing 6 > large enough, one shows that the truncated model is a good approximation 
of the original model, in the sense that the probability that they will differ by x goes exponentially 
fast to zero in x (Compare with (2.30) and (2.34) in [18j). 

Lemma 1.2.6 Let b = 6/a and r > Ee"^Vlog2. Then 



X . 



[Lr{uj) — Lr{uj) > x) < 2exp ^ ^ 
Proof: Fix 6 > and a positive integer r > 1 {N = r^). Notice that 



N 



< Lrico) -Lr{u)) < ^ZiI{Zi > blogr} . 



1=1 



By Markov's inequality, 



N 



N 



On the other hand, 



g|Zi/{fZi>f logr} 



qJ, aZ-\ 

= gf Zi/{e"^i>rTe^} 

< 1 + e^^^Ile"^^ > r-e—} 

„aZ\ 



(1.2.13) 
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and hence, 



Now, 



N 



^ZiI{Zi > 61ogr} > X < e~t^ 



log 



1 + 



ab 
T 2 



1 + 



ah 

r 2 



N 



< r 



log 1^1 + 



Ee 



aZi 



ah 

r 2 



a6 
7" 2 

4-a6 



< log 2 . 



if we take b = 6/a and r > Ee'^^^/log2. Together with ()1.2.13p . this proves Lemma 11.2.61 



□ 



Lemma 1.2.7 Under (ll.2.4p . there exist constants bi > 0, such that for all r >bo 



u 



P (|L (0, (r, r)) - EL (0, (r, r)) | > n) < 6i exp -b2j. , (1.2.14) 

V {^ogr)^yrJ 

for n G (0, 63 r'^/^ log r] . 

Proof: We have checked all the conditions of Lemma [1.2. 31 applied to the truncated process, where 
we take xq = 4(6 log r)^r, cq = 61ogr, ci = 1 and C2 = 1/4(6 log r)^ (this follows from Lemma ll.2.51 
and (|1.2.12p ). So there exist 03,04 > such that 

P (M, - M„ > «) < {2 + 1} exp (-c.^j^^^) . 

for ah u < 2(6 log r)r^/^. 

Using Lemma 11.2.61 and the fact that Lr > Zj.,we can see that there exists M > 0, such that 
for all r > 1, |E(L^) - E{Lr)\ < M. Therefore, for u> M, 

P(|L,. -E(L,.)I > 2n) < ¥{\Lr-'E{Lr)\>2u-M) 

< ¥{\Lr -'E{Lr)\ > U) 

+ F{\Lr -Lr\>U-M) 

= P {Mn -Mo>u) 

+ ¥{Lr-Lr >U-M). 

Again using Lemma ll.2.6( we can choose 60 = Ee"^i/log2, 61 > large enough and 62,63 > 
small enough such that (|1.2.14p holds not only for u G [2M, 63 (log r)r^/^], but also for < u < 2M. 

□ 
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Lemma 1.2.8 There exists a constant 64 > such that 

7r - 64(log^ r)y/r < EL (0, (r, r)) < -/r . (1.2.15) 

Proof: We note that the right hand side of ()1.2.15p follows from the definition of 7. To prove 
that the left hand side of p.2.8p also holds, we parallel the arguments developed by Howard and 
Newmann in [17J. 

We start by noting that it is enough to prove (I1.2.8P for integer values of r. Denote by Hr the 
set of points (x, t) such that x,t > and x + t = r. Then 

L(0, (2r, 2r)) < max L(0,x)+ max L(x, (r, r)), 

and hence, by symmetry with respect to H^, 

EL (0, (2r, 2r)) < 2E max L(0, x) . 

Define, for k € {0, 1, . . . , 2r}, x^. = (k, 2r — k), and for k G {1, . . . , 2r}, = (/c — 1, 2r — k). For 
any x G Hr, there exists a A; G {1, ... , 2r} such that < x. We have 

L (0, x) < max {L(0, Xfc_i) , L(0, x^)) + L(zfc, Zfc + (1, 1)}. 

This implies that 

EL(0, (2r,2r)) < 2E max L(0,Xfc) + 2E max L(zfc, z^ + (1, 1)). 

0<k<2r l<fc<2r 

The second term on the righthand side is bounded by the expectation of the maximum of the 
total weight in the 2r squares [z^, z^ + (1, 1)], which is clearly bounded by clog r, for some constant 
c > (here we use (jl.2.4p ). So we get 

EL(0, (2r,2r)) < 2E max L(0, x^) + clog r. 

By dm}, 

max EL(0,x) = max EL (O, (V2r — xv/x, \/2r — xJHc)) 
= EL(0,(r,r)) , 

and thus 

EL(0, (2r,2r)) < 2E max L(0, Xfc) + clog r 

< 2E max {L(0, Xfc) - EL(0, x^)} 

0<fc<2r 

+ 2 max EL(0,Xfc) + clogr 

0<fc<2r 

< 2E max {L(0, x^) - EL(0, x^)} 

0<k<2r 

+ 2EL(0, (r,r)) + clogr . 

(1.2.16) 
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Let 

Mr = max {L(0,Xfc) - EL(0,Xfc)} . 

0<k<2r 

Define for a large constant C > 0, tlie event 

2r 

A=f] {L(0, Xfc) - EL(0, Xfc) < C7(log2 r)^}. 

fc=0 

Tlien 

< C(log2r)V^U+L(0, (2r,2r))lAc. 

Tlierefore, 

EM^ < C(log2 r) + a/E [L(0, (2r, 2r))2] .p(y4^). 

We crudely bound L(0, (2r, 2r)) by the total weight in the square [0, (2r, 2r)], and see that there 
exists a constant c' > such that 

E [L(0,(2r,2r))2] < c'r^ 
We can use Lemma 11.2.71 to conclude that 

2r 

P(^'^) < J^P(L(0,Xfc)-EL(0,Xfc) >C7(log2r)Vr) 

fc=0 

< (2r + l)6iexp(-62Clogr). 
By increasing C, this shows that there exists c" > such that for all r > 1, 

EMr < c"{log^r)^/^. 
Together with ()1.2.16p . this proves that there exists 64 > such that for all r > 1 

EL (0, (2r, 2r)) - biilog^ r)^ < 2EL (0, (r, r)) . (1.2.17) 
By Lemma 4.2 in [17J, (|1.2.17p implies that the left hand side of (|1.2.15p is true. 

□ 

Proof of Theorem ll.2.2t The results of Lemma 11.2.71 and Lemma 11.2.81 now easily combine to 
Theorem 

□ 

1.3 Geodesies and a- Rays 

The second subject of interest to us are the geodesies. The existence of semi-infinite geodesies (or 
rays) for percolation like models has already been extensively study by Newman and coauthors 
[2.1]. They developed a general approach, based on Theorem 11.2.21 and on the curvature of the limit 
shape, that leads us to what they called the (5-straightness of geodesies. This property is the key 
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for proving the existence of rays. 

For each p G and 9 G (0, -7r/4), let Co(p, 6) denote the cone through the axis from to p and 
of angle 6. Let i?o"*(p) be the set of points q > p such that p G zi7(0, q). For fixed d G (0, 1), we 
say that the geodesies starting at in Co((l, 1),^) are (5-straight if there exist constants M, c > 
such that for all p G Co((l,l),6') and ||p|| > M, 

i2r(p)eCo(p,c|pr^). (1.3.1) 

Our argument on how to control the fluctuations of the geodesies will very closely follow the 
proofs given in [23] for the classical Hammersley process. In the classical case, the control of 
L(0, (r, r)) around its asymptotic value is stronger than our Theorem II. 2. 21 but our result is strong 
enough to extend the method to the more general Hammersley process. For details of the proof, 
we refer to |23j . 

For each L > 0, let 5Cyl(p, L) denote the side-edge of the truncated cylinder of width L, that 
is composed of points q G with q > p and |q| < 2|p|, and such that the Euclidean distance 
between q and the line through and p equals L. Assume that p = {x,t) G Co((l,l),^), that 
q e R7\p) and that q G 5Cyl(p, |p|^-''). Then 

L(0,q) =L(0,p) + L(p,q) 

or, equivalently, 

/(q) - /(q - P) - /(P) = A(0, p) + A(p, q) - A(0, q) , 
where /(p) = f{x, t) = ^\fxi is the shape function and 

A(p, q) = -L(p, q) - /(q - p) . 

Since qG 9Cyl(p,|p|i-^), 

/(q)-/(q-p)-/(p)>c|p|l-2^ 

for a finite constant c > 0, depending on 6* G (0,7r/4) (here we use the curvature of the shape 
function; see Lemma 2.1 in [23]). Notice that if 5 G (0, 1/4) then 1 - 2(5 G (1/2, 1) and so 

|p|l-25 » |p|l/2log|p| _ 

Hence, by Theorem 11.2.21 for b G (0, 1/4), we must have that if q G 9Cyl(p, |p|^^''), then with very 
high probability q This can be formalized to prove the following lemma: 

Lemma 1.3.1 Fix b G (0,1/4) and 6 G (0,7r/4). For each p = {x,t) G Co((l,l),6') and q G 
9Cyl(p, IpI^^''), /etG5(p, q) he the event that that there exists p' G p+[0, 1]^ such thatp' G uj{0,q). 
Then there exist constants K,Ci > such that, if \p\ > cq then 

P(G5(p,q)) <Cle-^2|Pl^ 

We extend this Lemma to hold uniformly for q and p in a fixed-size finite box, then use the 
boxes around q to cover the side-edge of the cylinder to get: 
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Lemma 1.3.2 Fix 5 G (0,1/4) and 9 e (0,7r/4). For each p = {x,t) G Co((l, 1), 0), let Gs{p) be 
the event that there exists p' G p + [0,1]^ and q G 9Cyl(p, |p|"'^^'^) such that p' G tu(0,q). T/ien 
there exist constants K,Ci> such that, if \p\ > C3 i/ien 

P(G5(p)) <C4e-^5|p|\ 

Now we can show 5-straightness by "gluing" together these cyhnders: if a geodesic starts close 
to p, with high probability it will exit the bottom edge of the cylinder Cyl(p, |p|^~''). Then we 
cover this bottom edge with boxes p2 + [0,1]^, where |p2| > 2|p|, and for each of these p2 we 
consider the cylinder Cyl(p2, |P2r~'')) and so on. With Borel-Cantelli we can make the probability 
that a geodesic starting close to p will ever leave through the outer edges of the boundary cylinders 
very small. The cylinders at the next step of the procedure have a slightly different angle than the 
cylinders in the previous step, but the changes in these angles are bounded by a geometric series, 
which means that all cylinders are contained in a cone starting at 0, of angle IpI^*^- This is basically 
the same argument used for the proof of Lemma 2.4 in [23] in the classical set-up. This leads us 
to: 

Lemma 1.3.3 Fix 5 G (0, 1/4) and 9 G (0, 7r/4). There exist constants n,Ci > such that for all 
p G Co((l,l),^) with IpI > cg, we have 

P I I U ^S"*(P ) I C Co(p, IpI-*^) ) > 1 - CTe-'^^IPl^ . 
V \p'ep+[o,ip / / 

Furthermore, with probability one, there exists M > such that for all p G Co((l,l),0) with 

IpI > M, 

i?r(p)cCo(p,|p|-^). 

An up-right (down- left, resp.) semi-infinite path starting at x is the lowest continuous increasing 
path through an ordered sequence (pn)n>o in I^^) with po = x, p„ G P and p„, < pm (pn > Pm, 
resp.) whenever n < m. We call (p„)„>o a semi-infinite geodesic if w{pn,Pm) C (Pn)n>o (every 
part of the path is a geodesic). Finally, for each angle a G (0, 7r/2) U {tt,3tt/2), we will call a 
semi-infinte geodesic (pn)n>o an a-ray if 

lim = a := (cos a, sin a) . (1.3.2) 

n— >co Ip^I 



Theorem 1.3.4 With probability one, every semi-infinite geodesic is an a-ray for some a G 
(0,-7r/2) U (7r,37r/2), and for every x G and a G (0, 7r/2) U {tt,3tt/2) there exists at least one 
a-ray starting at x. 

Uniqueness and coalescence of a-rays do not depend upon (5-straightness. The proof of these 
can be done by using a method introduced in [21], that would work in a wide context. In [23], 
Wiithrich applied this method to the classical Hammersley mode to get uniqueness and coales- 
cence for fixed directions. Here we state without proof the analogous result for the Hammersley 
model with random weights. The reader can convince her- or himself of the validity of the the- 
orem by checking that the proof given by Wiithrich can be adapted mutatis mutandis to our set-up. 
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Before we state the next theorem, we shah define what we mean by convergence of paths: we 
say that a sequence of paths tu" converges to w, and denote lim„_^oo = if for ah bounded 
subsets B cM? there exists uq such that w"^ Ci B = w Ci B for all n > uq. 

Theorem 1.3.5 For fixed a G (0, 7r/2) U (7r,37r/2), with probability one, for each x G there 
exists a unique a-ray starting at x, which we denote by ci7q,(x), and if 

nm - — - = a , 

n— >oo ||z„|| 

then 

lim tu(zo,z„) = CTa(zo) . 

Furthermore, for any x, y G there exists c„(x, y) G M? such that tUQ,(x) and zUaiy) coalesce at 
Ca(x, y). 

1.4 Busemann Functions 

In the middle of the fifties Busemann [4j introduced a collection of functions to study geometrical 
aspects of metric spaces. These functions are induced by a metric d, and by a collection of rays 
(semi-infinite geodesies) as follows: the Busemann function with respect to a ray {w{r))r>o, 

is the limit of 

•) := d{w{r),w{0)) - d{w{r), •) 

as r goes to infinity. Along a ray w the metric d becomes additive. By using the triangle in- 
equality, this implies that the defining sequence is nondecreasing and bounded from above, and so 
it always converges. Using analogous considerations, one can construct Busemann functions over 
spaces equipped with a superadditve "metric" L (one needs the reversed triangle inequality). 

Using the concept of a-rays, we will study the function Ba{x, y), which is defined by taking the 
first coalescence point c = c(a, x, y) between the a-ray that starts from x and the one that starts 
from y (remember that these two rays coalesce), and setting 

S,(x,y)=L(c,y)-L(c,x). (1.4.1) 

Note that if we take a different coalescence point c', then c > c' and they both lie on a geodesic. 
Since L is additive on a geodesic, we get 

L(c',x) = L(c',c) -FL(c,x), 

which shows that the definition of i?a(x,y) does not depend on the choice of the coalescence point. 

Let (z„)n>i be any unbounded decreasing sequence that follows direction (cos a, sin a), and 
let c(z„,x,y) denote the most up-right coalescence point between ■a7(z„,x) and ti7(z„,y). By 
Theorem 11.3.51 with probability one, there exists uq > such that 

Vn>no c(z„,x,y) = c(Q;,x,y) 



and so 



Vn>no L(z„,y) - L(z„,x) = B„(x,y) 



(1.4.2) 
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Therefore, in geometrical terms, i?a(x, •) can be seen as the Busemann function along the ray ^^(x). 

Some properties of B^ are summarized in the following propositions: 
Proposition 1.4.1 The distribution of the function B^ is translation invariant: 

S„(- + z,- + z) = S„(-,-). (1.4.3) 

The Busemann function is anti- symmetric and additive: 

Ba(x,y) = -Ba(y,x) and Ba(x, z) = ^^(x, y) + 5„(y, z) . (1.4.4) 
Fix X, y G and p, q G such that p, q > 0. The function 

Ah^ S„(x + Ap,y + Aq) (1.4.5) 

is cddldg in A G M. 

Proof: The translation invariance of the underlying compound two-dimensional Poisson process, 

{wq : q G (P + z)} = {wp : p G P} , 

implies (jl.4.3p . Anti-symmetry follows directly from the definition of the Busemann function. Now, 
by taking a coalescence point c between ^^(x), ^^(y) and ^^(z) we have that 

i?a(x, z) = L(c, z) — L(c, x) 

= -^^(c,z) - L(c,y) + L(c,y) - L(c,x) 
= -Ba(y,z) -hSa(x,y) , 

which clearly shows additivity, and finishes the proof of (|1.4.4p . 

By additivity and anti-symmetry, to prove cadlag, we can restrict our attention to y = (0,0), 
A = and q = (1,0) so that e is varying on the horizontal direction close to the origin. For 
different values of y the argument is similar. Choose c as the coalescing point of vOa{0) and ^^(q). 
Furthermore, define y <0 such that (0,y) is the crossing point of n7o(q) with the y-axis. Then for 
every AG [0, 1], we have that 

5«(0,Aq) =L(c,Aq) -L(c,0). 

This is because the a-rays starting at Aq are wedged in between ^^(O) and ^^(q). Furthermore, 
for A small enough, there will be no Poisson points in the rectangle [0, A] x [y, 0], which means that 
for those A, L(c, Aq) = L(c, 0), and therefore 

5,(0, Aq) =0. 

This proves right continuity. The existence of the left limit follows from monotonicity. 

□ 

We also have the symmetries: 
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Proposition 1.4.2 For any {x,t) G R^, 

B5n/2^a{0, (x, t)) ^ B^{0, {t, x)) (1.4.6) 

and 

Ba{0, {x, t)) = B^{0, {rx, t/r)) , (1.4.7) 

where tana = r~^tana. 

Proof: Define S as the reflection in the diagonal x = t. Then S(P) has the same distribution as 
P, and X < y < S{y). This shows that for all x, y G M^, 

L(x,y)^L(5(x),S(y)). 

Note that the lowest geodesic in the reflected case is not necessarily the reflection of the original 
lowest geodesic. However, we can now use (jl.4.2p to prove (|1.4.6p . The scaling invariance p.4.7p 
follows directly from symmetry (ll.2.3p . 

□ 

We now prove integrability of Ba- Before that we require the following lemma: 

Lemma 1.4.3 Let y = 0, a = bn/A, and x = (—1,-1). Let z be the intersection point between 
^577/4(0) o.'iT'd the one dimensional boundary of {p : p < x}. Then 

EL(z,0) < 00. 

Proof: Assume that c > is a small constant and that |z| < cr. Then 

L(z,0) < min {L ((-1, -cr), 0) ,L((-cr,-l), 0)} < P(cr) 

where P{r) denotes the total weight of the compound Poisson process in the set [—1, 0] x [—cr, 0] U 
[-cr,0] X [-1,0]. Hence 

P (L(z, 0) > r) < P (|z| > cr) + P (P(cr) > r) . 

By choosing c small enough, one can make P (P(cr) > r) integrable over r > 0. It is therefore 
enough to prove that P(|z| > cr) is integrable over r > 0. We first consider the case where z < 
(— 1, — cr); the case z < {—cr, —1) can be handled in exactly the same way. 

By Lemma 11.3.31 ((5-straightness of geodesies): fix (5 G (0,1/4) and 6 G (0,7r/4), then for all 
p G Co((l,l),0) with IpI > cg, we have 

U (P) I C Co(p, |p|-^) I > 1 - C7e-^«IPl" . (1.4.^ 

p'ep+[o,i]2 

Let 

b := (tan(237r/16))"\ b' := (tan(2l7r/16))"^ and p„ := {-bn, -n) . 

Fix no > big enough so that Co(pn, IPnl"*^) does not intersect {{—b't, —t) : t > 0} for all n > uq 
(recall that 6 G (0, 1/4) is fixed). Since the asymptotic direction of ^5,^/4(0) equals (—1, —1) it will 
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eventually be to the left of {{—b't,—t) : t > 0}. Therefore, if z < (— 1,— cr) and cr > no, there 
must be an n > cr and G (p^ + [0, 1]^) n ^5^/4(0) such that 

i2r(Pn)5Z:Co(p„,|p„r^). 

Using ()1.4.8p . we see that 



■(z < (-1,-cr)) < cje 



■C8|p„|'= 



n>cr 



for cr > max{no,C6}, which shows integr ability. 

□ 



Proposition 1.4.4 //x < y and x 7^ y then 

-Ba(x,y)>0 and <EBa{^,y) < 00. (1.4.9) 

Proof: Since L(c, y) > L(c,x) whenever c < x < y we have that Ba{x,y) > whenever x < y. 
Note that 

B,(x,x + (1,1)) > L(x,x + (1,1)). 

Clearly, EL(x, x + (1, 1)) > (there is a positive probability that a Poisson point will fall in between 
X and X + (1, 1)). Also, 

EB„(x, X + (1, 1)) = E5^(x, X + (1, 0)) + ES„(x, X + (0, 1)) . 

Property ()1.4.6p then shows that 

E5„(x,x + (1,0)) =ES^(x,x+ (0,1)) >0. 

Properties (I1.4.3P and (|1.4.4p show that for h > 0, 

EBaix, ih,0)) = /iEB„(x,x + (1,0)) > 0, 

and the same for the vertical direction. 

To prove that its expected value is finite, without loss of generality, assume that y = 0, that 
a = 57r/4, and that x = (—1,-1). Let z be the intersection point between ^5^/4(0) and the one 
dimensional boundary of { p : p < x}. By taking a coalescence point c < z < x, we have that 

B5,/4(x,0) = L(c,0) -L(c,x) 

= L(c,z) + L(z,0) - L(c,x) 

= L(z,0)-{L(c,x)-L(c,z)} 

< L(z,0). 

By Lemma ll.4.31 this proves integrability. 

□ 

The most important aspect of Ba is a Markovian structure described below: 
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Proposition 1.4.5 For all s <t and x E M we have 

5„((0, s), (x, t)) = sup {B„((0, s), (z, s)) + L((z, s), (x, t))] . 

z<x 

Proof: Without loss of generality we can take s = (and therefore t > 0). Define = Zaix, t) E M 
as the crossing-point of the a-ray starting at (x, t) with the x-axis. Clearly, Za < x and 

B^iO,{x,t)) = S„(0,(Z„,0)) + 5„((Z„,0),(x,0) 
= Ba{O,{Za,0)) + L{{Za,0),{x,t)). 

The last equality follows from the fact that (x,t) and {Za,0) are lying on an a-ray. This means 
that it is enough to prove that for all z < x, 

BaiO, {z, 0)) + L{{z, 0), (x, t)) < Ba{0, (Za, 0)) + LiiZa, 0), (x, t)) . (1.4.10) 

Suppose p is a coalescence point of the a-rays starting at 0, {x,t) and (z,0). Then 

i?«(O,(z,0)) =L(p,(z,0))-L(p,O) 

and 

B^{0, (Z,, 0)) = L(p, (Z«, 0)) - L(p, 0) . 
Furthermore, since p, (Zq,0) and (x,t) are elements of Wa{x,t), we know that 

L(p,(z,0)) +L((z,0),(x,i)) < ^(P,(a:,t)) 

= L(p,(Z«,0))+L((Z,,0),(x,t)). 

From this, ()1.4.10p easily follows. 

□ 

For t E M, define the positive measures on M, using Proposition 11.4.11 by 

ui{{x,y]) := B^iix, t), (0,0)) -Ba{{y,t), {0,0)) 

= B^{{x,t),{y,t)), (1.4.11) 

where x < y. Proposition 11.4.5) then shows that the process (f^)tgiR is a Markov process in the 
space of positive and locally finite measures defined on M: the future evolution of z^^ depends on the 
Poisson process in the upper-half plane M x (t,oo) and on the present value of z^* , not on the past 
of the process (which is of course independent of the Poisson process in M x (t, oo)). By translation 
invariance (Proposition 11.4. ip the distribution of does not depend on t: 

^i = <- (1-4.12) 

Thus, this distribution is an equilibrium (or time invariant) measure for the underlying Markov 
process. In the next section we will describe the generator of this Markov process, which will be 
an extension of the classical Hammersley interacting particle process. 



Chapter 2 

The Hammersley Interacting Fluid 
System 



2.1 The Graphical Construction 

It is well known that the classical Hammersley model, where all weights are 1, described in p], has a 
representation as an interacting particle system. The Hammersley Last Passage Percolation model 
with random weights has a similar description, although a better name might be an interacting 
fluid system. 

We start by restricting the compound Poisson process {ujp : p G P} to R x Then we choose 
a positive, locally finite measure v defined on M. Usually, these measures will be purely atomic, 
but this is not necessary. To each measure v we associate a non-decreasing process z^(-) defined by 



Note that i'(-) is a cadlag function. Although the details are a bit cumbersome, all the results we 
will show can be extended, mutatis mutandis, to the case where = — oo for x < 0, which would 
correspond to a non-locally finite measure with an infinite fluid density to the left of 0. This is a 
quite natural starting condition, but we will not use it explicitly. 

The Hammersley interacting fluid system (M*)j>o is a stochastic process with values in the 
space of positive, locally finite measures on M. Its evolution is defined as follows: if there is a 
Poisson point with weight a; at a point (xo,t), then 



Here, M*~ is the "mass distribution" of the fluid at time t if the Poisson point at (xo,t) would 
be removed. To the left of xq the measure does not change. In words, the Poisson point at (xo,t) 
moves a total mass u to the left, to the point xq, taking the mass from the first available fluid to 




for X > 
for X < 0. 



(2.1.1) 



M*({xo}) = Af*-({xo}) + u; 



and for x > xq 



M*((xo,x]) = (Af*-((xo,x])-^)+. 



(2.1.2) 
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the right of xq. See Figure \2A\ for a visuaUzation in case of atomic measures, of the process inside 
a space-time box. 

It is not true that the evolution M* is well defined for all measures (e.g. if we start with a 
finite number of particles to the left of 0, every particle would be pulled instantaneously to — oo). 
We will follow the Aldous and Diaconis [1] graphical representation in the Last Passage Percolation 
model (compare to the result in the classical case, found in their paper): 

Proposition 2.1.1 Let N he the set of all positive, locally finite measures v such that 

liminf^>0. (2.1.3) 

y->-oo y 

For each v G N , the process defined by 

L^{x,t) ■.= sup{iy{z) + L{{z,0),{x,t))} (x G M, t > 0) (2.1.4) 

z<x 

is well defined and the measure 

Miiix, y]) := LAy, t) - L,(x, t) , (2.1.5) 
evolves according to the Hammersley interacting fluid system. 

Proof: It follows immediately from the definition that L^[x,t) is increasing in x and t, even if Ly 
would not be finite everywhere. This implies that if we can prove that Ly is finite on for example 
Z X Z+ with probability one, then almost surely, L^ is finite everywhere. Therefore, we need only 
prove for any fixed point {x,t) that Ly{x,t) is finite with probability one. 



Use (|1.2.3p to see that for any e > and r = y^\x — z\t we have that 

P (L {{z, 0), (x, t)) > -ez) = P (L ((0, 0), r, r)) > -ez) . 
Using Theorem 11.2.21 we get 

¥{L{{z,0),{x,t)) > -ez) < P(L((0,0),(r,r)) >e(-x + r2/t)) 

for some positive constants c and c', and for large enough r. Borel-Cantelli then shows that 

L{{z,0),{x,t)) 
—z 

as z ->■ -oo. Now ([2X3]) gives us ([2X1]) . 



^0, 



We can show that x i— )• Ly{{z, 0), (x, t)) is cadlag using the fact that for y > x, 
L((z, 0), (y, t)) < L((z, 0), (x, t)) + L((x+, 0), {y, t)) 

and 

limL((x+,0),(y,t)) =0, 

yix 
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(0,t) 



{x,t) 



t/2 
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(0,0) 
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Figure 2.1: In this picture, restricted to [0,x] x [0,t], the measure consists of three atoms of 
weight 5, 3 and 7. The compound Poisson process consists of two points of weight 4 and 7. The 
measure MJ consists of three atoms of weight 1, 4 and 6, while at time t, it consists of one atom 
with weight 7. A total weight of 4 + 6 has left the box while a total weight of 2 has entered. 
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where x+ means that you are not ahowed to use a possible Poisson point directly above x. There- 
fore, M* is indeed a locally finite measure on M. 

To see that Af* follows the Hammersley interacting fluid dynamics that we have just defined, 
suppose that there is a Poisson point in {xo,t) with weight uj. For x < xq, this Poisson point has 
no effect, so M* = M*" on ( — oo,a;o). Clearly, 

Ml{{xo}) = Mt-{{xo}) + oj. 

If x > Xq, then the longest path to {x,t) that attains the supremum in (j2.1.4|) can either use the 
weight in {xo,t), which would give 

Ljy(x, t) = Liy{xQ,t) = L^{xo, t-) + UJ , 

or it could ignore the weight in (xo,t), which would give 

Ly{x,t) = Lu{x,t-) . 

This proves that 

M*((xo,x]) = L^(x,t) - L^(xo,t) 

= max{L,^(xo,t-) + oj , Ly{x,t-)} - L„{xo,t-) - uj 
= (M*-((xo,x])-^)^. 

□ 



2.2 The Generator 

Recall that we have defined the Hammersley interacting fluid system M* as a Markov process on 
the space N (Proposition l2.l!T]) . To compute the generator G of M*, we endow the space N with 
the weak topology of convergence of continuous functions with compact support, so t'n — >• iff for 
all (j) e Cc(M), we have 

/ (f){x)vn{dx) — / (j){x)v{dx). 
Jr Jr 

Now assume that we have a continuous bounded function / : A/" i— )• M that only depend 
on "cylinders" , so on the configuration of the fluid in a compact interval. Consider the transition 
operator 

Ptf{u) := E {f{M') \M^ = u) . 



Proposition 2.2.1 Let IZz^w represent the operator acting on M that moves a total mass uj to the 
left, to point z, taking mass from the first available fluid to the right of z. Let f : M M be a 
continuous bounded cylinder function. Then 

t^o t 
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where 



oo fOO 



oo JO 



[f{RzA^)) - f{y)] dF{oj) dz . 



Proof: Suppose / only depends on the configuration v restricted to the interval [a, b]. We wish to 
determine Ptf{i') for small values of t. Choose large with N > —a, and define Xi{t) > X2{t) > . . . 
as the x-coordinates of the rightmost Poisson point in the strip (— oo, — A^] x [0,t], the second 
rightmost point, etc. Define uJi{t) as the weight of the Poisson point with x-coordinate Xi{t). 
The configuration of M* at [a, b] will not depend on the compound Poisson process in the strip 
(—00, a — N] X [0, t], if we have the event 



E, 



^ojj{t) < u{[Xi{t),a)) Vi > 1 



Define the following random walk in 



Snit) 



p2 . 



-Xi{t) + a- N,Y,^iii) 



i=l 



The steps to the right are iid with distribution Exp(t), and the steps up are iid with the distribution 
F of the weights. We define the corresponding step- function for x E [0, 00): 



i=l 



Clearly, is a decreasing family (as t — t- 0) of increasing functions. In fact, for fixed t, x^*^ 

describes a compound Poisson process, with intensity t and jump distribution F. Define 

ip{x) = i^{[a — N — X, a)). 

Since G M, we know that v satisfies (|2.1.3|) . so there exists e > and A'^o > 1 such that for all 
X > and ah N > Nq, 

'ij){x) > ex + eN. 

The event Et occurs when — "^l^) foi^ all x > 0. Since the distribution of the weights 

F has an exponential moment, we know that there exist constants ci,C2 > (depending only on e 
and Nq) such that for all t small enough and > A'o, 



P(V X > : x^^\x) <ex + eN)>l- cite 



-C2N 



This follows from the classical result where the weights (i.e., the jumps up) have an exponential 
distribution. We conclude that 

F{Ef) < cite-"^^. 



Now choose N^ = — log(t). Also, define the event 



At = {the strip [a — Nt, b] x [0, t] contains or 1 Poisson point}. 
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Then, given that M° = I/, we get 

E(/(M*)-/H) = E(li5,nA,(/(M*)-/H)) 

+ K{lE^uAi{fiM')-fiu))) 
= E(/(M*) - f{u) I Et n At){l - ¥{Ef U A?)) 
+ E{lE^uA^^{f{M')-f{u))) 
t{b-a + Nt) 
l + t{b-a + Nt) 

roo -| pb 

-/(i/)] dzdF{uj) 

+ o(t) . 

The o(t)-term follows from the fact that f{E^) = o{t), P(A^) = o{t) and / is bounded. The 
factor in front of the integrals equals the probability of having 1 Poisson point in the strip, given 
that the number of points equals or 1. Of course the x-coordinate of this one point will have a 
uniform distribution on [a — Nf, b]. The proposition now follows when we divide by t, take the limit 
t — 7- 0, and note that Poisson points to the right of b have no influence on f{M^^^), so 

fb poo 

/ f{RzA^)) - f{y) dz = / f{R,M) - f{u) dz. 

Ja-Nt Ja-Nt 

□ 

Assume that we have a probability measure \E' defined on M and consider G as a realization 
of this probability measure. We say that is time invariant for the Hammersley interacting fluid 
process (in law) if 

Ml = M^ = v for all t > . 

In this case, we also say that the underlying probability measure on N is an equilibrium measure. 
In particular, a probability measure ^' on is an equilibrium measure if it satisfies 

j Gf{v)d-^{v) = 0. (2.2.1) 

It was this interplay between the longest path description and the equilibrium interacting parti- 
cle system that proved very fruitful in the results for the classical Hammersley process in [6l [7| . We 
will attempt the same in the interacting fluid system, but since the equilibrium solution to (I2.2.ip 
is not explicitly known (in most cases), we needed to develop new tools and ideas, which in fact 
also had interesting applications for the classical case. 

2.3 Exit Points 

We have two remarks concerning the interactive fluid process. The first remark concerns the flux of 
the fluid system through the t-axis. We define the flux measure v* on [0, oo) such that for t > s > 0, 



i/*((s,t]) = L,(0,t)-L,(0,s). 



(2.3.1) 
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In Figure [2J] we can see that z/*((0,t]) = 10, whereas i^*((0,t/2]) = 4. For x,t >0, we also have 
the equahty 

L^{x,t) = max < sup {iy{z) + L{{z,0), {x,t))} , 

^ 0<z<x 

sup {u*{z)+L{{0,z),{x,t))}} . 

0<z<t ^ 

(2.3.2) 

The way to prove ()2.3.2p leads us to our second remark: the supremum in the definition of 
Lu{x,t) is actually attained (it is therefore a maximum). This is the statement of the proposition 
below. 

Proposition 2.3.1 With probability 1, the set {z < x : L^{x,t) = v{z) + L((z, 0), (x, t))} is 
non-empty, for all {x,t) G M x [0,cxd). Furthermore, if we define 

Z^{x,t) = sup{z < X : L^{x,t) = iy{z) + L{{z,0),{x,t))}, (2.3.3) 

then 

L^{x, t) = u{Z^{x, t)) + L((Z^(x, t), 0), (x, t)). 
Also, the following local comparison property holds for x < y and t > 0: if Z^{y,t) < then 

L^{y, t) - L^{x, t) < L(0, (y, t)) - L(0, (x, t)) ; 

If Z^{x,t) > then 

L(0, {y, t)) - L(0, (x, t)) < L^{y, t) - L^(x, t) . 

Proof: The first part of this proposition relies on the fact that on a compact interval, the sum of a 
non-decreasing right-continuous function (z ^ z^(^)) and a non-increasing left-continuous function 
{z I—)- L{{z,0), (x,t))) attains its maximum, and the set of maxima is closed. 

Now we have to show that with probability 1, for all {y,s) G M x [0,oo), the supremum over 
z < y can actually be restricted to a compact set. For a single (x, t) (and therefore for a countable 
set of (x, t)'s), this has been done in the proof of Proposition 12. 1 . 1] using (j2.1.3p and Theorem ll.2.11 
for L. Then we can conclude the desired compactness property for all {y, s) by using the inequality 

u{z) + L{{z, 0), {y, s)) < u{Z,{x, t)) + L{{Z,{x, t),0), {y, s)) 

for all z < Zy{x, t), y > x and s < t. 

This last inequality can be seen as follows: suppose there exists z < Zy[x,t) such that 

u{z) + L{{z, 0), (y, s)) > u{Z,{x, t)) + L{{Z,{x, t),0), {y, s)). 

Then the geodesic zu{{Z^{x,t),0), (x,t)) must intersect the geodesic w{{z,0), (y,s)) in some point 
c. It follows that 

L{{z, 0), {y, s)) = L{{z, 0), c) + L(c, {y, s)) , 

and 

L((Z^(x, t), 0), (x, t)) = L{{Z^{x, t),0),c) + L(c, (x, t)) . 
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Furthermore, we have that 

i^iz) + L{{z, 0), {y, s)) > u{Z^{x, t)) + L{{Z^{x, t), 0), c) + L(c, (y, s)). 

Combining this gives 

v{z) + L{{z,{)),{x,t)) > u{z) + L{{z,0),c)+L{c,{x,t)) 

> u{Z^{x,t)) + L{{Z^{x,t), 0),c) 

+ L(c, {x,t)) 

= u{Z^{x,t)) + L{{Z^{x,t),0),{x,t)) , 
which contradicts the definition of Ziy{x,t). 

Now we prove the second statement. If Zi,{y,t) < then there exists Zy < such that 

Lu{y,t) = L{{zy,0), {y,t)) + ly(zy). 

Define z as the (an) intersection between the geodesic from {zy,0) to {y,t) and the geodesic from 
to {x,t). Then we can see that 

Lu{x,t) > L{{zy,0),{x,t)) + iy{zy) 

> L{{zy,0),z) + L{z,{x,t)) + u{zy). 

This imphes that 

L^{y,t) - L^{x,t) < L{{zy,0),{y,t))-L{izy,0),z)-L{z,{x,t)) 

= L{z,{y,t))-L{z,{x,t)) 

= L(z,(y,t))-L(0,(x,t)) + L(0,z) 

< L{0,{y,t))-L{0,{x,t)). 

The proof of the third statement fohows the same line. 

□ 

Note that we have also proved the following statement: 

Vy>xVO<s<t : Z^{y, s) > Z^{x,t). (2.3.4) 

Equation (j2.3.2p now can be seen as follows: define i^- as the restriction of i' to (— oo,0] (so has 
no mass on the positive x-axis). Then for x,t>0, there exists Z <0, such that 

L,_{x,t) = i^iZ) + LiiZ,0),{x,t)) . 

The finite geodesic w{{Z,0), (x,t)) crosses the positive t-axis in some point, call it (0,s). It is not 
hard to see that 

u*{s) = u{Z) + L{{Z,0),{0,s)) 

and 

L,_ {x, t) = u*{s) + L((0, s), {x, t)) > u*{s) + L((0, s), (x, t)) 
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for all < s < t. This follows from similar arguments as at the end of the proof of Proposition 12. 3.11 
Finally, we remark that 

Ly{x,t) = max < sup{i/_(2:) + L((z, 0), (x, t))} , 

2<0 

sup{z.(z) + L((z,0),(x,t))}| , 

2>0 J 

from which Equation (|2.3.2p follows. 

2.4 Law of Large Numbers for Exit Points 

Define the following crossing points: 

{Za{x,t),Q) is the crossing point of Wa{x,t) and M x {0}. (2-4.1) 
This was already used in the proof of Proposition 11.4.51 Analogously, 

{0,Z*{x,t)) is the crossing point of uJaixjt) and {0} x M. (2.4.2) 
Recall p. 4. lip and let fa := u^. We have the following proposition. 
Proposition 2.4.1 For all x & M and t > 0, we have 

Za{x,t) = Ziy^[x,t). 

Proof: From ()1.4.10p we immediately get that Z^^{x,t) > Za{x,t). Now suppose Z^^{x,t) > 
Za{x,t). Since 

Ba{0,{x,t)) = L^^{x,t) 

= Ua{Z,^ {x, t)) + L((Z,„ (x, t), 0), (x, t)) 

= B^{0, {Z,^ {x, t), 0)) + L{{Z,^ (x, t), 0), (x, t)) , 

we conclude that 

L((Z,„ (x, t), 0), (x, t)) = Ba{{Z,^ (x, t), 0), (x, t)). 

Now define c as the coalescing point of the a-ray starting at {Z^^ (x, t), 0) and the a-ray starting 
at (x,t) (which goes through (Zq,(x, t), 0)). Since the two points are on one a-ray, then 

L(c,(x,t)) = Ba{c,{x,t)) 

= S„(c,(Z,Jx,t),0)) +5«((Z,Jx,t),0),(x,i)) 
= L(c, {Z,^ (x, t), 0)) + L{{Z,^ (x, t), 0), (x, t)) , 

which would imply that there exists a longest path from c to (x, t) which is strictly below the 
original a-ray, contradicting the uniqueness of the lowest finite geodesic. 

□ 
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Since the a-ray WaiO) has a.s. the asymptotic angle a, by translation invariance, we can easily 
see that 

lim P(|Z„(t, (tana)i)| > et) = lim ¥(\Z*{t, {tan a)t)\ > et) = . (2.4.3) 

However, we can also control for more general v. Compare the following theorem to Lemma 
3.3 in tl4J. 

Theorem 2.4.2 Suppose v S M . Assume that 

liminf^^^Ml > 2^/^ and limsup ^^^^^^ < JVt^. (2.4.4) 

Let /i G M. Then, with probability one, 



Zyit + h, (tan a)t) 
lim = , 

t— >oo t 



Furthermore, define 



Z^^h{t,{i8.na)t) = argsup{z^(z) + L(-t„ + (z,0),(/i,0))} . (2.4.5) 

z<t 

where := (t, (tana)t). Here we take the right-most maximum. The idea is that we move the 
origin to —t^ and look at the exit point for {h,0). Then, with probability one, 

Zi, h{t, {tan a)t) 
lim ^ = . 

t— >CXD t 

Proof: We start with the proof for Z,^. Using the transformation (|1.2.3p . we can assume, without 
loss of generality, that a = bir/A. 

The proof of this lemma is based on the following elementary estimate for the shape function 
f{x, t) = jy/xi: for s > 

' }s-^^ ifs<8t, 
f{t + s,t)-f{t,t)<{ i{s>St. (2-4.6) 

Fix e > and suppose that Z^^h{t) := Z^{t + h,t) < —et — 1. Define k = — \Zi,{t)'] and n = \ t + h\. 
It follows that 

u{-k + 1) + L{{-k,0),{n + l,n + 1)) - L((0, 0), (n, n)) >0. 

Choose r/ > small enough (we will see how small). For t big enough, we know from (j2.4.4p 
that ^ 

^{—k + 1) < + 

Here we use that k > et. This implies that for t big enough, 

L{{-k, 0), (n + 1, n + 1)) - L((0, 0), (n, n)) > ^jk - r]k . (2.4.7) 
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A straightforward application of Theorem 11.2.11 and Borel-Cantelli shows that for t big enough 
(and therefore for n big enough), we will have that 

|L((0, 0), (n, n)) - /(n, n)| < r/n. 



Below we will work out a more complicated application of Theorem 11.2.11 We now have for t big 
enough 

L[{—k, 0), (n + 1, n + 1)) — f{n, n) > -7/c — 7]k — rjn . 

For t big enough, we will have that k > ^en (since h is fixed). Now we subtract f{n + k, n) — f{n, n) 
from both sides of the inequality and use ()2.4.6p to see that 



is greater or equal to 



L((-A;,0),(n + l,n + 1)) - f{n + k,n) 



— ryfc — r/n if ^en < < 8n , 



■^=k — r]k — r]n if /c > 8n . 



which is greater or equal to 



£ 7 
128 



n — 9rjn li k < 8n , 



-^k - 2r]k if A; > 8n . 

k V 8 

Therefore, if the set {t > : Zy{t,t) < —et — 1} is unbounded, it follows that for some small 
rj > 0, the event that 



|L((-A;,0),(n + l,n + l)) -f{n + k,n)\ > < 



' rjn if /c < 8n , 
^ rjk if k > 8n . 



happen infinitely often for n > 1 and A; > 1. Using Borel-Cantelli, this will have zero probability if 
for all 7] > 0, 



00 8n 



^ J]] P {\L{{-k, 0), (n + 1, n + 1)) - f{n + k,n)\> r]n) < +00 (2.4.^ 



n=l k=0 



and 



^ ^ P {\L{{-k, 0), (n + 1, n + 1)) - /(n + k,n)\> r]k) < +00. (2.4.9) 

n=l k=8n 

Theorem 11.2.11 gives us some control on the fluctuations of L about its asymptotic shape. Note 
that for < /c < 8n, we have 

P(|L((-/c,0),(n + l,n + 1)) - f{n + k,n)\ > rjn) = 

P(|L(0,r(l,l)) - f{n + k,n)\ > r]n) . 
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; n 



for r = Y^(n + 1 + fe)(n + 1). 

Let u = v?!"^ . If we choose n large enough, we can make sure that u < r]n/2 and 

|/(n + k,n) — jr\ < r/n/2. 

Also, for n large enough, we have that u S [ca-^log^r, C4r^/^logr] (see Theorem ll.2.ip . This 
implies, using Theorem ll.2.1l that there exist q > such that for n large enough and < A; < 8n, 

P(|L(0,r(l,l)) -/(n + A;,n)| > r/n) < F{\L{0, {r,r)) - jr] > u) 

< ci exp (^—C2n^^^/ log 

This clearly proves (I2.4.8p . 

For k > 8n, we define u = /c^/^, and in a similar way we find that 

P(|L(0,r(l,l)) -f{n + k,n)\ > rjn) < ci exp (^-C2k^/^ / logk'^ 
This proves (l2X9]l . 

The proof that the set {t > : Zy(t + h,t) > et + 1} is bounded with probability 1 is actually 
easier, since we will have that < A: < n. Therefore, we only need to use the following bound on 
the shape function f{x,t) = 'jy/xi: for < s < t 

fit-s,t)<l,s-Yj- 
The remainder of the argument is similar to the previous case. 

Now we consider Zi,^h. Since for all t we have 

Z^,h{t,t) = Z^{t + h,t) 

(but not as processes!), the Borel-Cantelli type arguments for Z,^^h hold in this case as well, mutatis 
mutandis. 

□ 



Chapter 3 



Busemann Functions and Equilibrium 
Measures 

3.1 Busemann-Equilibrium Measures 

Fix a G (7r,37r/2), recall (jl.4.1ip . and define the measure 

:=z.O((x,y]) =5„((x,0),(y,0)) for x<yGR. 

Of course, as an immediate consequence of Proposition 11.4.51 for all G M x [0, oo), 

L,Jx,t) = 5,((0,0),(x,t)). (3.1.1) 

Now we are able to make the first step to characterize the equilibrium measures for the Hammersley 
Interacting Fluid process: 

Proposition 3.1.1 Fix a G (7r,37r/2) and denote by the probability measure on N induced by 
the random measure Ua- Then ^'^ is an equilibrium measure for the Hammersley Interacting Fluid 
Process. 

Proof: By (|1.4.12p and (fSXT]) . we have the following: 

MlJ{x,y]) = L,^{y,t) - L,^{x,t) 

= i?„((0,0),(y,i)) -i?„((0,0),(x,t)) 
= '^i{(.x,y]) 
= i'a{{x,y]), 

which shows that is indeed an equilibrium measure for the Hammersley Interacting Fluid process 
M*. 

□ 

We will show that the measure I'a has the following mixing property, usually called strong 
mixing in dynamical systems. We consider the cj-algebra J- = (T{z^Q,((a, 6]) : a < b £ M} on the 
sample space il., defined by the compound Poisson process restricted to M x M_. We can define the 
translation Tt as an J^-measurable map from Q to 0, simply by translating all Poisson points by 
the vector (i, 0). 
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Proposition 3.1.2 The random measure has stationary and integrahle increments. Further- 
more, it satisfies 

\/A,BeF: \\m ¥{Ar\T^'^{B)) =¥'{A)'^{B). (3.1.2) 
In particular, this implies that is ergodic. 

Proof: Stationarity and integrability follow from Proposition 11.4.11 and Proposition 11.4.41 From 
translation invariance and a standard approximation of sets in it is enough to prove (j3.1.2p for 
&\\ A,B eFh:=(y{vo,{{a,h]) : a<be[0,h]}. 

Consider the geodesies 

w{—ta, 0) and w{—ta, {h, 0)) . 

where ta ■= {—t,—t{tana)). Almost surely, these paths will converge to roa(O) and zi7q,((/i, 0)), 
respectively, on any finite box. This means, that if we define for a,b G [0, h] 

4*)((a,6]) =L(-t«,(6,0)) -L(-t„(a,0)), 

then for t big enough, we have = fa|[o,?t]- 

Clearly, T^f^{B) is independent of iJa\ since they depend on the Poisson process to the left 
respectively to the right of the line {—t} x M. Define the event 

Ct = {\fs>t: =i.„|[o,,,]} 

and denote A^^^ the counterpart of the event A in 

-^f :=cT{i^i*)((a,6]) : a<be[0,h]}. 

To see what is meant by j4(*\ we define the index-set 

I = {iya{{a,b]) : a<6e [0,/i]}, 

and 

/W={i.W((a,6]) : a<be[0,h]}. 

There is a canonical bijection i : I ^ I^^\ Define B as the product cr-algebra on M^, and likewise 
B^^h Extend the canonical map i such that i:W -i- M^^'' . Define the map : ^ by 

u t-^ (^{uj) := {iya{ia,b]){uj) : a<b£[0,h]}, 
and likewise (/>t : — t- M^^'' . We know that 

Fh = r\B) and fP = ^^HB^'^) ■ 
This means that there exists U £ B, such that A = (p~^{U). We define 

A^'^ = ct>T\im . 
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Thus 

\¥{AnT-_^\{B))-nA)F{B)\ < \F{AnT-^\{B)nCt) 

-F{A)F{B)\ 
+ P(Cf) 

= \F{A^'^nTf_^\{B)nCt) 

-F{A)F{B)\ 
+ P(Cf) 

< |P(AWnr-^\(i?)) 

-F{A)F{B)\ 
+ 2P(C7f ) 
= F(S)|P(AW) -P(A)| 

+ 2P(C7f ) 

< 4P(C7f). 

The proposition now follows from the fact tha10 P(Cf ) — t- 0. 

□ 



3.2 Uniqueness of the Equilibrium Measure 

Now we state the most important result of this chapter. 

Theorem 3.2.1 If we start the Hammersley Interacting Fluid Process with then 

M^" = for all t>0. 
The process x i— )• I'aix) is stationary and ergodic and its intensity is given by 

7(F) 



Ez^^(l) = iL^Vtana. (3.2.1) 

Finally, consider a random v G M , which is time invariant, and which defines a stationary and 
ergodic process on M with Ez^(l) G (0,oo). Define a G (7r,37r/2) by 

/ 2 , . 
a = arctan — -— Ez^il 



7(F) 



Then v'= v^. 



Proof: The first statement is an immediate consequence of Proposition 13.1.11 The fact that 
X I— 7- fa(x) is stationary, ergodic and that Et'a(l) G (0,oo) follows directly from Proposition 13.1.21 



^We note that, in the classical model, we have independent increments even if the probability of the event CI does 
not decay to very fast. This indicates that, to show mixing by using these events may not be the best strategy. 
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Fix a G {tt,3tt/2) and set p = p{a) := Vtana. Since the model is invariant under the map 

{x,t) ^ {px,t/p), 

Eua{l) = Ez^5^/4(p) = Ez^57r/4(l)Vtana. 

By Proposition 12.4.11 
Now, for all t > 0, 

Lu,^/At^t) = (Z5,/4(f)) +L((Z5,/4(t),0),(i,t)) . 

If ^5^/4(0 > then 

< L,,^/M^t) - L(0, < z.5,/4 (^5^4(0) ■ 

On the other hand, if Z^^i/^{t) < then 

zL/,{t) ■.= z;^^^^{t,t) = z;^^,{t,t)>o. 

It follows from ([2XT]1 and (f3XT]) that 

i^*^ix) = B^{0,iO,x)). (3.2.2) 
From the additivity of the Busemann function, we know that 

which finally implies that, 

< L,^^^^{t,t)-L{0,{t,t)) 

< max|i/5^/4(Z5^/4(t)),zy;^/4(Z*^/4(t))| . 

Clearly, from p.4.6p (Proposition I1.4.T]) it follows that 

Let e > 0. The ergodicity of i^^jr/i (and of ^^5^/4) and the fact that Efa(l) € (0,oo) imply that 
there exists r/ > such that for all t big enough, 

P (max{z.5,/4(r?t) , ^^5*^/4(7/*)} > rf) < e/2. 

By (I2X3D . for any rj > 0, 

P (max {z5,/4(t), Z5*^/4(t)} > 7?t) < e/2 , 
for large enough t. Combining all this gives, for t big enough, 

P(L,3^/,(t,t)-L(0,(t,t)) >et)<e. 



3.2. UNIQUENESS OF THE EQUILIBRIUM MEASURE 



39 



Therefore, 

t 

By (jS.l.ip , and using Proposition 11.4.11 we get 



0. (3.2.3) 



= Ei?5,/4(0, {t, 0)) + Ei?5,/4((t, 0), (t, t)) 
= 2iEK,/4(l)). 

The ergodic theorem applied to ^5^/4(0, {t,t)) implies thalU, with probabiHty one, 

hm = 2Ez^5^/4 1 . 

t— >-oo t 

Combining this with Theorem 11.2.11 and ()3.2.3p , one gets p.2.1|) . 

Now we need to address the uniqueness of f^. Suppose £ M is ergodic and time invariant. 
Define Z(t) = Zi^(t, {ttana)) and Zfi{t) = Z^(t + h, (tana)t). Now define, as in Theorem 12.4.21 

Z{t) = argsup {u{z) + L{-ta + (z, 0), 0)} 

z<t 



and 



Z,,(t) = arg sup {u{z)+L{-t^ + iz,0),{h,0))} . 

z<t+h 



Here, we take the right-most location of the maximum. The intuition for Z{t) and Zh{t) is that we 
place again the origin at —to,, and look at the exit-point for the path that starts at —to,, picks up 
mass from v and then goes to 0, resp. {h,0). 



Clearly, we have 



{Z{t),Zn{t))^{Z{t),Zh{t)). 



Since v is ergodic, and by our choice of a, v satisfies ()2.4.4p . Theorem 12.4.21 then tells us that 

t-\Z{t),Zn{t))^Q. 

This means that the two geodesies 

vj{-ta + {Z{t),{)),Q) and tz7(-ta 0), (/i, 0)) 
will converge, in any , to the a-rays 

vOa{0) and g7q((/i, 0)) 

respectively (this follows from Theorem ll.3.5p . However, these two a-rays will coalesce, which 
means that with probability 1, there exists to > such that for all t > to, the two converging paths 
coalesce, which in turn implies that Z{t) = Zh{t) (because they are both the right-most point where 



^The same method to prove Proposition 13. 1.2 1 can be used to prove ergodicity of -85^/4. 
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the maximum takes place and, as soon as they coalesce, they get the same exit point). 
Now define 



and 



We also have that 



L{t) = sup{u{z) + L{{-t + z,- [tan a)t),0)} 

z<t 



Lh{t)= sup {z^(z) + L((-t + z, -(tana)t), (/i,0))} . 

z<t+h 



{L{t),Lh{t)) = {L^{t, {tan a)t),L^{t + h, {tan a)t)). 



Furthermore, if t > toi then 

Lh{t) - L{t) = S«(0, {h, 0)) = i/«((0, h]). 

This proves that 

^t'{tan")(*^*'* + ^]) = L^(t + /i,(tana)t) -L^(t,(tana)t) 
A i/a((0,/i]). 

Since v is time invariant and ergodic, we see that 

z.((o, h]) ^ M(';,,„),((t, t + h])'^ ^.^((0, h]). 

In principle, we need to show convergence for a finite number of /I's simultaneously, but it is 
not hard to see that the ideas we used can be extended to that case, at the cost of some notational 
burden. 

Note that we have proved that for any deterministic v satisfying (I2.4.4p . 

as a process in h. This shows that in a rarefaction fan, the fiuid process converges locally to the 
correct equilibrium process (local equilibrium). 

□ 



Corollary 3.2.2 

E«(l))- ^(^^ 



2-v/tan a 

In particular, for all a G (vr, 37r/2), we have 

EK(1)).EK(1)) = ^, 
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Proof: Remember that 

u*^{x) = Ba{0, (0,x)). 

For a = 57r/4, the result follows from Proposition 11.4.21 and Theorem 13.2.11 Now use the map 

(x,t) I— 7- {px,t/p) to see that 

= 1^5^/4 (^x/ V tan aj . 

□ 



3.2.1 The Classical Hammersley Process 

Theorem 3.2.3 Consider the classical Hammersley model and let ^ p denote a probability measure 
on J\f induced by a one- dimensional homogeneous Poisson process of intensity p > 0. Then is 
an equilibrium measure. 

The proof of this theorem will be an immediate consequence of a version of Burke's Theorem 
(Theorem 13.2.41 below). To formulate this theorem, we start by considering M* on a fixed interval 
[0, i?]. Also, we choose z/* as an independent Poisson process of intensity 1/p. The evolution of 
restricted to [0, R] only depends on v* , on the Poisson process in the strip [0, R] x [0, oo) and on the 
starting configuration v on [0, i?]. This description is called the Hammersley process with sources 
{^\[o,R]) s-iid sinks (i^*). In fact, (Af^|[o,_R])t>o is still a Markov process, since z^* is a Poisson process. 

When we look at the paths of the Hammersley particles induced by M*, which in fact correspond 
to the level sets of L^{x,t), they enter the box [0,R\ x [0,T] either at the bottom (as sources, i.e. 
atoms of u), or at the east, i.e. the right-side of the box. We will call this process of "entries" Eji. 
Furthermore, the Hammersley paths exit the box on the left-side of the box as sinks (atoms of ly*) 
or at the top (atoms of MJ). Finally, the paths have bottom-left corners at the Poisson points of 
P, and we call the process of upper-right corners P*. 

Theorem 3.2.4 Using the notations introduced above, the process Er is a Poisson process of 
intensity 1/p, the process MJ is a Poisson process of intensity p and the process P* is a (two- 
dimensional) Poisson process of intensity 1. Furthermore, Er, M'^ and P* are all independent. 

The idea of the proof is to show reversibility of the process M*. We introduce the space E as 
the state-space of M*, so 

where Eq = {0} and 

En = {ixi,X2,. .. ,Xn) : < Xi < 2:2 < • • • < a^n < -R}- 

We endow E with the usual topology, which makes it into a locally compact space. We can view 
the measure as a measure on E. One of the goals of the following calculations will be to prove 
that v is an equilibrium measure. 

We can also write down the generator of the Markov process M*|[o^ij]: if / G Co{E), then for 
X £ E ^ 

Gfix) = [ f{lZtx)dt + -f{Cx) -{- + R] fix) 
Jo P \P / 
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where C corresponds to an exit to the left and TZt corresponds to an insertion of a new Poisson 
point at t. So C : E ^ E where 

^ J (X2, . . . , x„) if X £ En {n> 2), 
\ ifxeEoUEi. 

And for < t < R, TZt : E ^ E where 
TZtX = 

Here we use the convention that xq = 0. 



{Xi, . . . ,Xi-i,t,Xi+i, ...,Xn) if Xi-i <t <Xi {x e En) , 
{Xl,. . . ,Xn,t) if Xn < t {x £ En) . 



Note that there is a difference with Proposition 12. 2. H because of the sinks. Viewing G as an 
operator from L^{p) to L°°(//), we can calculate the dual operator G* . This will be the generator 
of the time-reversed Markov process. We define TZ as an exit to the right and Cg as a new point at 
s such that the point directly to the left of s moves to the right, that is TZ : E ^ E where 

^ J {xi, . . . ,Xn-i) if X £ En {n> 2), 
\ if X G ^0 U £^1 , 

and, for < s < R, C.:^ : E ^ E where 
Define 



xi,...,Xi_i,s,Xj+i,...,x„) if Xj < s < Xj+i (x e En), 

S,Xi,...,Xn) if S < Xi {x e En). 



G*9iy) = £giCsy)ds + -^giTZy) " ( J + ^) diy) ■ 

Lemma 3.2.5 For all f,gG Cc(M), 

Gf{x)g{x)d^,ix)= [ f{y)G*g{y)d^p{y). 

E Je 

Proof: See the Appendix in [B]. 

□ 



Proof of Theorem 13.2.41 Lemma 13.2.51 shows that G* is the generator of the time reversed pro- 
cess. Notice that G* exactly corresponds to the generator of a Hammersley process with sinks 
and sources that moves to the right: put the sinks on the right-hand side as a Poisson process of 
intensity and let the particles jump to the right to a Poisson process of intensity 1 in the plane. 

We now have two conclusion. First of all, since G*\ = 0, ^ p is an equilibrium measure for the 
classical Hammersley process. Secondly, since the time reversed process has the same generator as 
a Hammersley process moving to the right, which we denote by M, we conclude that the upper-left 
corners of the Hammersley paths P* must be a Poisson process of intensity 1, since they correspond 
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to the points the M particles jump to. 

The process of entry-points Eji is shghtly more subtle, since not all entries can be retraced by 
looking at the process M; after all, if M* = 0, there might be entries caused directly by sinks. 
However, this holds in exactly the same way for the right moving Hammersley process M, and 
therefore we can conclude that the process Eji, which are the entry-points for the Hammersley 
process, corresponds to the exit-points of the time-reversed Hammersley process, which in turn 
correspond to the sinks of the process M, which by construction are a Poisson process of intensity 
1/p. Furthermore, we can conclude that all three processes MJ, P* and E^j are independent, since 
this is true by construction for the process M. 

□ 

Corollary 3.2.6 Consider the classical Hammersley model. Then ^{Si) = 2 and Va is a one- 
dimensional homogeneous Poisson process of intensity p = Vtan a. 

Proof: Together with Theorem 13. 2. 1^ Theorem 13.2.31 implies that 

V / 2 

i^p = i^a with a = a(p) = arctan p 

We have seen that, for p = 1, z^i(l) = z^i(l) = 1- Therefore, by Corollarv 13.2.21 

4 

and hence 7((5i) = 2. In particular, p{a) = Vtana. 

□ 




3.3 The Multi-Class Process 

For two positive measures and on R, we say that dominates v, notation v > v, whenever 
^ foi' all measurable / CM. 

Proposition 3.3.1 Suppose we have two measures v^v ^ M such that P > z^. Define the corre- 
sponding interacting fluid system as and , using the same weighted Poisson process (basic 
coupling). Then Mj^ > (as measures). If i^\(-oo,o) = ^l{-oo,o)) then M^{[0,x]) — Mf{[0,x]) is 
non-increasing in t for all x > 0. 

Proof: Fix an interval [—K, K] and a time t. There exists (a random) M > such that 
and Mf restricted to [—K, K] only depend on Poisson points in [— M, K] x [0, t] and on and P 
restricted to [—M,K] (it is not hard to see that we can take M = Zi^{—K,t), see (j2.3.4p ). This 
means that we are only dealing with a finite number of Poisson points, so if we can prove that the 
premise "M]' > for all s < t" implies that > M^, we will have proved the first statement, 
since it is obviously true for t = 0. 
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Suppose there exists a Poisson point at {xo,t) with weight cj for some xo G [—M,K], since 
otherwise the imphcation is immediate. We then know, using Proposition 12.1.1] and (|2.1.2p . that if 
Xo <x<y, 

Mriix,y]) = (Mr_((xo,y])-o.)+-(Mr_((xo,x])-u;)+ 
< (Mf_((xo,2/]) -a;)+ - (M^ ((xq, x]) - 
= Mr{ix,y]). 

The inequahty follows from the fact that if A > B and A > B > 0, then {A — u)^ — (B — uj)-^- < 
(A + A — w)+ — (B + B — a;)+. If x < xq < i/ or x < y < xq, the implication is straightforward, 
following a similar split up. 

The second statement follows from a similar reasoning: suppose there is a Poisson point at 
(xo,t) with weight uj. If x > xq > 0, 

Mf((xo,x]) -M,-((xo,x]) = (Mt((xo,x]) 

-(M,-_((xo,x])-u;)+ 
< Mf_((xo,x])-M,1((xo,x]). 

The inequality follows from the fact that 

{A - c)+ -{B- c)+ <A+-B+ 
whenever c > and A > B. Since 

Mf ([0,xo]) - Mn[0,xo]) = Mt([0,xo]) - Mr_([0,xo]) , 

"t ills sllO'V^^S t tlctt 

Mf([0,x]) -Mr([0,x]) < Mt([0,x]) -Mr_([0,x]). 

Now suppose Xo < and x > 0. Note that under the condition on P, we have that for all s > 
and all e > 0, Lp{—e, s) = L^{—e, s), so 

Mf ([0, x]) - M^{[0, x]) = Lp(x, s) - L^(x, s). 

When Lp{x,t) does not use the weight at (xo,t), we know that Lp{x,t) = Lp{x,t—) and that 
Ljy(x, t) > Ly{x, t—), which implies the desired result. If Lp(x, t) does use the weight at (xq, t), then 
it is not hard to see that L^{x,t) will also use the weight at (xo,t) (the longest path corresponding 
to i> is always to the right of the path corresponding to v), which means that only the mass on the 
X-axis strictly to the left of is used, and therefore 

Mf([0,x]) = Mn[0,x]). 

Finally, when xq = x, we get that 

Mf([0,x])=Mt([0,x])+a; 

and 

Mr([0,x])=Mr_([0,x])+a;, 
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and when xq > x, we have 



Mf([0,x]) 



Mt([0,x]) 



and 



Mr([o,x]) 



MU[0,x]). 



□ 



In other words, Proposition 13.3.11 tells us that the Hammersley Interacting Fluid System is 
monotone: if one starts the fluid process with the same Poisson weights (basic coupling) and with 
ordered initial configurations, then the order is preserved for all t > 0. This coupled process is 
called the Multi-Class Fluid System. It is just a convention to describe a coupled process with 
ordered initial configurations |13j . 

For any countable D C (7r,37r/2), one can a.s. construct simultaneously a collection of equilib- 
rium processes {^a ■ ol S Z)} by using the same Poisson weights on M x M_ and the Busemann 
functions B^. It turns out that this collection respects the order induced by the angles q G D. 
More precisely: 

Theorem 3.3.2 If a > a then ^ ^a- In particular, for any countable subset {ui : z G Z} C 
if one runs simultaneously (basic coupling) the Interacting Fluid Processes on M x M-|_ with initial 
measures {ua- : i S Z) then, whenever Oi > aj, M^°'' > M^"^ for all t>0, and 



Proof: Almost sure coalescence of a-rays for fixed a G {tt,3tt/2) implies almost sure coalescence 
for a given countable D C (7r,37r/2). This allows us to construct simultaneously the Busemann 
functions {Bq. : a G D}, as a function of the underlying compound Poisson process. Thus, time 
invariance of the Busemann multi-class measure (f^j : i E Z) follows from translation invariance 
of the compound Poisson process (as in the proof of Proposition II .4. 1]) . 

To see that it is indeed a multi-class measure, let m be the crossing point between zj7q((z,0)) 
and Wai{z' ,0)). Furthermore, denote c as the coalescence point of the two a-rays Wa{{z,0)) and 
Waiiz' ,0)), and denote c as the coalescence point of the two a-rays vJaiizjO)) and Wadz' ,0)). 



(M^"^ : i G Z) = (i^«, : i£Z) 



Then 



{[z,z']) -i^4[z,z']) = {l(c,(z',0))-L(c,(z,0))} 



- {l(c,(z',0))-L(c,(z,0))} 
= L{c,{z',0)) 



-{L(c,m) + L(m, (z',0))} 
+ L{c,{z,0)) 

-{L(c,m) + L(m, (z,0))} 

> 0. 



□ 
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This result is also new in the classical Hammersley Interacting System, where a different and 
explicit description of the multi-class equilibrium measure with a finite number of classes is given 
by Ferrari and Martin [T3j . 



Chapter 4 

Second-Class Particles 



4.1 Second-Class Particles and Exit Points 

Proposition 13 . 3 . ll can be used to define the notion of second-class particles. In the interacting fluid 
system we can define it analogously to the interacting particle case, with a slight adaptation due 
to the continuous weights. 



We start by changing v into by putting an extra weight e > in 0, so 

P([0,x]) = z^([0,x]) + e for x>0. 

With this new process, and using the same Poisson weights, we define M/'. Clearly, 

Mf([0,x]) <M,'^([0,x])+£. 

Now define the location of the second class particle Xy{t^ as 

X^(t) = inf{x > : Mf ([0,x]) = M/'([0,x]) + e] . 

By Proposition 13.3.1] is a non-decreasing function of t, meaning that the second class 

particle moves to the right. In fact, the extra mass e will spread out, and our definition coincides 
with the rightmost point of this spread-out mass. This is a natural choice, since we will show that 
it does not depend on the total mass e, while for example the leftmost point does depend on e. 



There is the following important connection between the longest path description and the second 
class particle. Let be the process defined by v^(x) = v{x) for x > 0, and by ^^{x) = — oo for 
X < 0. We also define the process v~ by setting v~ [x) = for x > 0, and v~[x) = u{x) for x < 0. 
Then 

sup{L((z,0),(x,t)) +z^(z) :0<z<x} ifx>0 
— oo if X < , 



Ly+{x,t) 
and 



L^- (x, t) = sup{L((2;, 0), (x, t)) + v{z) : z < and z < x} . 

Clearly, 

Ly{x,t) = max{Lj,+ (x,t),L,^-(x,t)} . 
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Now suppose X > 0. If L^+{x,t) > Lj^-(x,t), there exists a longest path that does not use any 
weight of I' on (— oo, 0). This means that if we add a weight e > in the origin, 

Lp{x, t) = Lp+ {x, t) = Liy{x, t) + e . 

Using Proposition 12. 1 . H we see that this means that 

Mf(x) = Mnx)+e 

so X^{t) < X. 

If on the other hand we start with Xu{t) < x, we conclude that 

Mf(x) = Mnx) + e 

using Proposition 13.3.11 and the fact that and are right-continuous. This in turn means 
that 

Lp{x,t) = Ly{x,t) +e, 
which is only possible if L^+{x, t) > L^- (x, t). 

We have shown that 

{X,{t) <x} = {L,+ (x, t) > L,- (x, t)} . (4.1.1) 
Note that this can be rewritten as 

{X^{t) <x} = {Z^{x, t) > 0}. (4.1.2) 

This means that the path of the second class particle corresponds to a competition interface, a 
fact well known for the totally asymmetric exclusion process |15j . This allows us to show that the 
second class particle satisfies a strong law whenever i/^ and have asymptotic intensities. The 
proof of this does not use a coupling of two invariant versions of the fluid process, as is usual in the 
interacting particle case, but it uses the longest path description in a direct way. 

We would like to point out that in our general set-up, with random weights on the Poisson 
points, we do not have an equivalent of Burke's Theorem. This means that the time-reversed 
process is not a Hammersley interacting fluid system. Therefore, the path of a second class particle 
in general does not coincide in law with a longest path in the interacting fluid system, in contrast 
to the classical case, where the statement is true [7]. However, we do have the following connection. 

Proposition 4.1.1 Assume that the distribution ofv is translation invariant. Then, for any t > 0, 
we have that 

X^{t) - X = -Zy{x,t). 

Proof: This follows almost immediately from ()4.1.2p . since that equality can be rewritten as 

{X^(t) -x<h] = {Z^{x + h,t) > 0}. 
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Now use translation invariance to see that 

Zv{x + h,t) = Zy{x, t) + h. 

□ 

When we consider all a-rays starting at the line M x {t} and we move from left to right, {Xy^ {t) , t) 
is the first point where the a-ray passes the origin. It is tempting to think that the a-ray starting 
at {Xy^{t),t) actually passes through the origin, but this is false in general. In fact, after time t, 
most a-rays will have coalesced with other rays, and the crossings with the rr-axis will be quite far 
apart; we would conjecture they are order t^/'^ apart. 



4.2 Law of Large Numbers for Second-Class Particles 

Proposition 14.1.11 allows us to use Theorem 11.2.11 and Theorem 12.4.21 to prove a strong law for 
the second class particle in the case of Va- However, we are able to prove a strong law even for 
deterministic initial conditions that satisfy a density property given in the following lemma, whose 
proof is very similar to Theorem I2.4.2[ 

Lemma 4.2.1 Suppose v £ N . Assume that 

hminf inf "^^^ " ^' " + > ^ Vti^ . (4.2.1) 

d-^oo ze[-Ad,Ad\ 2d 2 

Let e > 0. Then, with probability one, 

hminfMl±£M^>0. 

t— >oo t 



Now assume that 



u{{z- d,z + d]) 7 



limsup sup -^^^ ^ < T^Vtana. (4.2.2) 

Then, with probability one. 



ZJ (1 - e)t,( tan a)t) 
limsup — — < 0. 

i— >co t 

Remark: Note that in our density condition, we do not allow the midpoint of the interval to be 
much larger than d. The reason for this might be more clear if we think of a Poisson process: if we 
fix d, we can always find some z G M such that the interval [z — d, z + d] is empty! 

Proof of Lemma 14. 2. It The proof of this lemma relies on the proof of Theorem 12.4.21 As usual, 
we will assume without loss of generality that a = 57r/4. We will start with the first statement. 

We need to prove that for all > small enough, the set 

{t > : Zy{{l + 2e)t,t) < r]t} 
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is bounded with probability one. Choose t] < e/2. Define the translated measure iy2et such that for 
all b > a, 

i^2et{{a, b]) = u{{a + 2et, b + 2et]) . 
Using the definition in (j2.1.1[) . we see that for z G M, 

i^2et{z) = u{z + 2et) - u(2et) . 

Define k = - \Z^{{1 + 2e)t, t) - 2et] and n= [t\. lit is such that 

Z^((l + 2e)t,t) <r/t, 

then 

u{-k + 1 + 2et) + L{{-k + 2et, 0), ((1 + 2e)t, t)) 

- [u{2£t) + L((2ei, 0), ((1 + 2e)t, t))] > . 

This event has the same probability as the event 

i^2et{-k + 1) + L((-fc, 0), (t, t)) - L((0, 0), {t, t))>0. 

We know that k > et (since r] < e). We can use ()4.2.ip to see that for t big enough, 

-iy2et{-k + 1) = iy{{-k + 1 + 2et, 2et]) > ]pk - rjk . 

To see this, define half the length of the interval hy d = {k — l)/2, and the midpoint hj z = 
et — {k — l)/2. We see then see that \z\ < Ad. This gives us, using the fact that n < t < n + 1, 

L{{-k, 0), (n + 1, n + 1)) - L((0, 0), (n, n)) > ]pk - rjk . 

This is exactly ()2.4.7p in the proof of Theorem 12.4.21 ci^id we can follow that proof from this point. 
The second statement of the lemma follows similarly. 

□ 

With this lemma we can proof the following result: 
Theorem 4.2.2 Assume that 

lim — ^ = lim -—^ = - Vtana. (4.2.3) 

X— >oo X a;— >— oo x 2 

Then, with probability one, 

X,{t) 1 



lim 

t-!-oo t tan a 

Proof: As usual, we will assume without loss of generality that a = 57r/4. Suppose £ > and 
X^{t) <t — 2et. Define n = [t\. Then for t large enough, we have 

X^{n) < (1 - 2e)n + 1 + 2e < (1 - e)n . 
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By (f4T2]) . this implies that Zu{{l - e)n,n) > 0. 



If we can show that v satisfies (I4.2.2p . then we can use Lemma [4.2.1l to see that Z^{{l—e)n, n) > 
can happen only for finitely many n > 1, which gives 

liminf^^ > 1. 

t— >oo t 

Bounding the limit from above can be done using the analogous argument. 



To see that (I4.2.2P indeed holds, remark that for all 1 > > 0, there exists R > 1 such that for 
ah z> R, 



z.((0,z]) 1 



< r] and 



z.((-z,0]) 1 



< T]. 



Choose M > 1 such that R\) < M. Now choose d > MR/i]. U z > d + R or z < -d - R, 

we get that 

v{{z-d,z + d]) < d'y + {2\z\ + 2d)r] 
< dj + Wdr] . 

If \z + d\< R, 



u{{z-d,z + d]) < -{2d + R)'j + M + \z - d\T] 
< d-f + {5d + 2 + ^j)r, , 

and a similar bound holds when — d| < i?. This proves that satisfies (j4.2.2p . 



□ 



4.3 The Rarefaction Fan 



Assume now that we have a non-homogeneous initial measure u that satisfies the following rarefac- 
tion assumption: 



1- v^ix) ^ [y) 

au ■= limsup < limmf =: o^, 



x^oo X y^-oo y 

In this regime the behavior of the second class particle is different from the previous one [lOj . 



(4.3.1) 



Theorem 4.3.1 For each v € 



// ()iXT]l holds then 



7^ 7^ 



let 



a{v) := arctan(l/?;) G (0,7r/2) . 
> as t ^ oo , 
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where Vu is a random variable with support 



7^ 7^ 
46^' 4^ 



and 



'{V,<v)=F[ snpMz) - (z)} > snpMz) - u^^,) (z)} ] . (4.3.2) 

\2>0 2<0 w / 

for V G ( jrz , jT" ) , where Va o,nd u are independent random measures. 



Proof: Condition ()4.3.ip ensures that with probabiUty one, at least one of the two suprema is 

(2 2 
IP") 4^ 

that at the boundary, (j4.3.2p may fail!) 



For each t>, t > let 

a = = arctan(l/i;) G (0, 7r/2) , 
X = xt := t/ tana and '■= {x , x tan a) . 

(j4.1.ip can be rephrased to 

{X^{t) <x} = {L^+{x, t) - L(0, (x, t)) > L^- {x, t) - L(0, [x, t))} . 

Thus, 

{X^{t) <x} = {L^+M > L^-M} 

sup [l^{z) + L{{z,0),Xa)] 
0<z<x 



>sup [i/(z)+L((0,z),x„)]| 

2<0 J 

sup [z^(2;) + L((2;,0),Xq,) - L(0,Xa)] 

0<2<X 

sup [v{z) + L((0, z),Xa) - L{0, Xq)] i 



> 

2<0 



(4.3.3) 



Notice that, for any compact set K C 



(L(z,Xq,) - L(0,Xq,) ; z G if) ^ (5a(0,z) ; z G i^) as x oo . (4.3.4) 

This remark is the core of the proof. It follows from Theorem 1 1.3. 5 1 we can take n > 1 big enough, 
so that K C [— n, n] x [—n,n]. The a-rays starting at {—n,n) and {n,—n) will coalesce at some 
point c G M^. Furthermore, the longest paths to x^, starting at (— n, n) and (n, —n), will converge 
to the respective a-rays in a bigger bounded square, containing c. From that time on, for all z £ K, 
we will have 

L(z,Xq,) - L(0,Xq) = Sq,(0,z). 

Notice that, since Ua is a function of the underlying two dimensional compound Poisson random 
set P restricted to the upper half plane, it is independent of i^. Together with (j4.3.3p and Lemma 
14.3.21 below, this proves the theorem for Z^/tan a G {au,bu)- 
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Now, if bu < 2 ^ tan a then a.s. 

supi/"(y) -v^[y) = <x). 
y<o 

By the same reasoning as in the proof of Lemma 14.3.21 

'D 

sup {1^(2;) + L((2, 0), Xq.) — L(0, Xq,)} — > 00 , as x — )• cxd . 

z<0 

Hence, both functions in (j4.3.2p are zero. The case where Ojy > ^-v/taiTa follows from a similar 
argument. 

□ 



Lemma 4.3.2 // ^Vtan a G {uu, bu) then, as x ^ 00, 



sup {v{z) + L((z,0),Xq) - L(0,Xq,)} — )■ sup{z^(z) - Va{z)} , 

0<z<a; i;>0 

and 

sup {1/(2;) + L((2;,0),Xa) - L(0,x«)} — ;> sup {1/(2;) - Va{z)} . 

z<0 z<0 



Proof: We will prove that < ^\/tana implies the first statement. The proof of the second 
statement follows exactly the same reasoning. 

Pick a' £ (0, 7r/2) and a' > such that 

< a = tan a' < ^ V tan a . 

Let Zai{x) be the crossing point between Wa'{0) and the line {{y,t) : y £ M} (recall that t = 
x/tana). Let be the event that Za'{x) is to the right of x^,. By Proposition 12. 3. T] (local 
comparison), we have that, under the event E^, 

L(0,Xa) - L((0,2;),Xa) > i^a'iz) 

(translate the origin to Xq and reflect both coordinates), and hence 

u{z) - [L(0,x„) - L((0,2;),x„)] < u{z) - v^'iz) . 

From Theorem 12.4.21 (recall that a' < a) , we conclude that 

lim W{E^) = 1. 

By our choice of a', a' > a^, and thus there exists M > (random) such that 

< Snp{l'{z) — Ua'iz)} = sup {I'iz) — Ua'iz)} < OO . 
z>0 0<z<M 
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This means that for each e > 0, we can find R > such that for all x large enough, the probability 
of the event that 

sup {u{z) + L{{z,0),Xa) - L{0,Xa)} 
0<z<x 

> sup {1/(2;) +L((z,0),x„) - L(0,Xa)} 

0<z<R 

is smaller than e. Now we can use (j4.3.4p to finish the proof of the first statement of Lemma I4.3.2I 

□ 

In the classical Hammersley process is a Poisson process of intensity 

1 



Iv 



which allows us to compute the law of the asymptotic speed in some cases. 
4.3.1 Poisson Initial Configuration 

Choose A > and fi > X, and assume that {u~^{y) ; y > 0) and {u^ [y) ; y < 0) are independent Pois- 
son counting processes of intensity A and ;U, respectively. We will only consider the case p G (A,^), 
the other cases are trivial. 

Define two asymmetric simple random walks and W~ , with 

F(W+(n + 1) - W+(n) = +1) = p+ := 

X + p 

and 

F{W- {n + 1) - W- (n) = +1) = p- :-- ^ 



p + p 

Since i^"*", ly^, v~ and v~ are independent Poisson counting process, it is not hard to see that 

sup{z^^(2;) — vti'^)} = supVF~'~(n) , 

2>0 n>0 



and that 



sup{i^ {z) — v p {z)} = swpW (n) 

2<0 n>0 



Furthermore, it is well known (and easy to see) that 



S' := snvW\n) ~ Geo(r^) , 

n>0 



where 

r' = 

1 -p' 
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since < 0.5 (for i = +,—). We find 

oo 

F{S^>S-) = ^P(5+ > fc)F(5- = A;) 

oo 

= 5^(r+)^(r-)^(l-r_) 
fc=0 

1 — r~ fJ- — P 
1 — r+r~ /i — A 

Therefore, by Theorem 14.3.11 (and taking = l/^/v), this proves that 



FiV > v) 



' if u > 1/A2 

1 -A 



if <t; < l/A^ 



1 if V < 

This agrees with the results found in [5l [12] . 
4.3.2 Periodic Initial Configuration 

Choose A > and /i > A, and assume that {i^^{y) y > 0) and (z^~(y) ; y < 0) are deterministic 
periodic configurations of intensity A and /i, i.e., i/"*" is concentrated on {k/X : > 1} while z/~ is 
concentrated on {A;/// : k < —!}■ Suppose X < p < fi. Theorem 14. 3 . 1 1 then tells us that we should 
consider 

P ( sup{i/(z) - Up{z)} > sup{u{z) - Up{z)} ] , 
\z>0 z<0 J 

where Vp is a Poisson process on R of intensity p. 
Define 

5"*" = sup{i^(z) — Vp{zy\ and S~ = sup{i^(z) — Uplz)}. 

2>0 z<0 

Clearly, for z > 0, we have ^{z) = [Xz\ , where [x\ denotes the integer part of x. Therefore, 

5+ = sup{-jyp{z) + [Xz]} 

z>0 

= [sup{- p{z) + Xz}\ 

z>0 
z>0 

The infimum of a Poisson process minus a linear function is studied in [22] . Theorem 3 in [22] 
entails that for k > 

P(5+ > A;) = ¥{\M{iyp{z) - Xz}] < -k) 

2>0 

= P(inf {i/,(z) - Az} < -A;) 
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where p+ = p+{X, p) G (0, 1) is the positive solution of 

For z < we have z^(z) = — [^|2;|J. So when we switch to positive z, we get 

= sup{up{z) - [fiz\} 

z>0 

= \snp{iyp{z) - fiz}]. 

z>0 

Now we can use Equation (7) of [22j: for A; > 

F{S- <k) = P(sup{z/p(z) - fiz} < k) 

z>0 

i=0 

In particular, ¥{S^ = 0) = 1 — p/fi. Therefore, 

oo 

P(5+>5-) = Y.¥{S- <k)p^{l-p+f 

k=0 

= p+C^-p/p) 

oo k 

X 

k=0 i=0 



We were not able to simplify this formula significantly. For A = 1 and /i = 2, Figure [^7\] gives 
the graph of P{S'^ > as a function of p G (1, 2). We compare it with the probability 

P- P 

pL — \ 

which we would get if we would take a Poisson process of intensity A for x > 0, and of intensity 
for X < 0. 

An interesting limit is A — )• oo. This means that = 0, and we get 

P(5+ > S-) = I- p/p. 

This is exactly the same probability as when we would have a Poisson process of intensity to the 
left of 0, and no particles to the right of 0! 



THE RAREFACTION FAN 




1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 
Figure 4.1: Picture of P{S'^ > S~), together with uniform distribution. 
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Chapter 5 

Fluctuation Results 



5.1 Some Elementary Identities 

We have seen that there is only one family of ergodic equilibrium measures for the Hammersley 
interacting fluid system. Let us denote it by {i^a : A > 0}, where 

A := Ei/a(1) = ^Vtana (5.1.1) 

(recall (fTO) ). Let 

(j)--(SLy = ^_ and V:=\- (5-1.2) 

^ V2A/ tana ^ 2A ^ ^ 

Thus, the speed (p corresponds to the characteristic angle a associated to the equilibrium measure 

ux. For simple notation, put Lx := Ly^, 

vx{x) := Ba{0, (x, 0)) , uiix) := 5«((0, t), {x, t)) , 

u*xit) := B^{0, (0, t)) , zv*'"(t) := 5,((x, 0), (x, t)) . 

Proposition 5.1.1 For given A > and x,t>0 

Yar Lx{x,t) = Varz^^(t) - Yar i^x{x) + 2Cov {i^xix), Lx{x,t)) . (5.1.3) 

In particular, 

Yar Lx{(t>t,t) = 2Cov {ux{(t>t), Lx{cl)t,t)) . (5.1.4) 
Proof: By Proposition 11.4.1] (additivitv). 

Lx{x,t) = i^x{x) + i^l'^'it) . 

Thus, 

Yar Lx{x,t) = Var z^a(3;) + Var i/^'^(t) + 2Cov (1/^(2^), '^A'^(i)) 
= Var ux{x) + Var i^^'^(t) 

+2C0V {iyxix),Lx{x,t) - i^xix)) 
= Varz^^'^(t) — Yav i^xix) + 2Cov (vxix), Lx{x,t)) 
= Yavul{t) -Yaruxix) + 2C0V {iyx{x),Lx{x,t)) . 
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(In the last step we also used translation invariance.) 

By Proposition 11.4.2] (rescaling). with 

a := 57r/4 and r := </)~^/^ = V tan a , 

we have that 

ux{(l)t) = B^{0, {4>t, 0)) = 55^/4(0, {(l>'/h, 0)) 

and 

Again, by Proposition 11.4.2] (reflection). 

i?5./4(O,(</)^/'t,O))=i?5./4(O,0^/'t) 

which shows that 
In particular, 

Vari^A(0i) = Varz^^(t) 

and hence, 

Yai Lx{(pt,t) = 2Coviuxi4't),Lx{(pt,t)) . 

□ 

Recalling that 

Lxix,t) := sup{ux{z) + L((z,0), {x,t))} , 

one may guess that the correlation between Lx and ux is proportional to the amount collected by 
Lx in the positive x-axis. Which is to say that 

Gov (uxix), Lxix, t)) ~ EZ+{x, t) (5.1.5) 

Since (p = (tana)~^, one may expect that 

where ^ G (0, 1) is the critical exponent for transversal fluctuations of a-rays. This indicates that 
the variance of Lx along the characteristic speed (f) has sub-linear growth. 

Another remarkable identity is given below: 
Proposition 5.1.2 For x G M and t > let 

bx,t ■■= X- <j)t. 

Then 

E{{Lx{x, t) - [yx{hx,t) + i^t]f ) = Var(LA {4>t, t) ) . (5.1.6) 
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Proof: By Proposition I1.4.1] (additivity). 

Lx{x, t) = iyxibx,t) + Ba{{bx,t, 0), (x, t)) . 
On the other hand (translation invariance), 

5a((&x,t,0),(x,t)) =S„(0,(x-6^,f,i)) =S„(0, t)) = LA((/>t,t). 

Since 

we get that 

La(x, t) - [vx{bx,t) + i^t] = Lx{(l)t, t)) - , 

which shows (|5.1.6p . 

□ 

To illustrate the importance of (j5.1.6p . let us assume that (jS.l.Sp holds and that zt is a path 
with speed a ^ (j), and let bt = zt — 4>t. By ()5.1.6|) . we get that 

Lx{zt,t)-[ux{bt) + 4^t] 

lim p = U , 

t^oo 

in the \? sense. Therefore, if one has Gaussian fluctuations for the equilibrium measures (which is 
a plausible assumption) then Lx will have the same behavior along any speed a ^ (p. This will be 
made precise in the classical Hammersley model. 

5.2 The Exit Point Formula 

Theorem 5.2.1 Consider the classical Hammersley model. Then 

Cov {vx{x),Lx{x,t)) = m[Z+{x,t)) . 

In particular, 

Var La(x, t) = ^ - Ax + 2AE (^^(x, t)) . 

Proof: For each e > we couple u'^ with i^A+e and let the negative part stay the same. Denote 
this new measure by vx,e- Let 

/(6):=E(L,, ^(x,t)) , 
and, for each n > let (disregarding the dependence on x and t for a while) 

a„ := ¥.{Ly^ Xx,t) \ uxA^) = • 

Note that does, in fact, not depend on e since we are condition a Poisson process on the interval 
[0,x]. Then 

^^^^^gx(A + er^_.(A-|..)^^ 

n=0 



62 



CHAPTERS. FLUCTUATION RESULTS 



and 



ae A ^-^ n\ ^-^ n\ 

n=0 n=0 



E(L,,z.a)-xE(L,J , 



which shows that 



A-/(0) = E(L,,z.a)-AxE(L,, 



= E(L,,z.a)-E(z.a)E(L,J 

= Cov(zyA,-^-A) • 



Now we are going to calculate this derivative in a different way. Note that the probability of 
finding more than one Poisson point of v\ ,^ — v\ lying within [0,x] is o(e). If we have exactly one 
Poisson point z G [0,x], then this point will contribute to Li,^^ if and only if z G [0, Zx{x,t)]. So 

/(e) = E(L,,,(x,t)) 

= E{Lxix,t))+e [ F{Zxix,t) > z)dz + o{e) 







Thus, 



= /(0) + e/ F{Zx{x,t)> z)dz + o{e). 
Jo 

^/(O) = 1^ {Zx{x, t)>z)dz = E {Z+{x, t)) , 



which shows the first part of the proposition. The second part follows by combining the first part 
together with Proposition 15.1.11 



□ 



Corollary 5.2.2 Let a G (7r,37r/2) and for x,t>0 define 

B^{x,t) :=B^{0, {x,t)). 

Then 

t 



V&iBa{x,t) = —f= — (V tan a)x + 2(V tan a)EZ(^(x, i) 
Vtana 

In particular, 

VarBa {t, (tana)t) = 2(\/tan a)EZ+(t, (tana)*) . 
Proof: It follows from Theorem 15.2.11 and Corollary 13.2.61 



□ 
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5.3 Gaussian Fluctuations 

To illustrate the importance of (jS.l.Gp . let us restrict ourselves to the classical Hammersley model. 
Thus, 

1 2 

4> '■= J2 ^^'^ ^ A ' 



Theorem 5.3.1 Consider the classical Hammersley model. Given a deterministic path {zt)t>o, ls.t 

bt = zt- X-h . 

ThcTi 

^^ E{{L,{z„t)-[.,{h) + 2X-H]}') 

t— s>oo t 

Proof: Notice that 

< t-^z+{x-^t,t) < x-^ . 



Together with the Dominated Convergence Theorem, Theorem 12.4.21 implies that 

E(Z+{X-\t)) 
lim ^ ^ ^ = . 

t— >CXD t 

By Proposition 15.1.2) and Theorem 15. 2. 1^ this proves Theorem 15.3.11 



□ 



Corollary 5.3.2 Consider the classical Hammersley model. Assume that {zt)t>o is a deterministic 
path such that 

lim — = a . 

t—>-oo t 

Then 

hm i = a := \aX — — . (5.3.1) 

2 



Furthermore, if a ^ X then 



Lx{zut) - [Xzt + X-H] 



lim F ' ' ' '—f i <u] = F{N < u) , (5.3.2) 



where N is a standard Gaussian random variable. 



Proof: Theorem 15.3.11 shows that 

Lx{zt,t)- [iyx{bt) + 2X-h] _ 

in the sense. Since i^x is a one-dimensional Poisson process of intensity A, this implies (|5.3.ip 
and (fO:2]l . 
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□ 



It is also interesting to remark that Proposition 15.1.21 implies that 

Li{zt,t) - [ui{bt) + 2t] V Li{t,t)-2t 



(5.3.3) 



Baik and Rains [3] proved that the right hand side of (j5.3.3p converges in distribution to Fq, the 
zero-mean Tracy- Widon distribution. Therefore, the left hand side of (j5.3.3p has the same limit 
behavior. 

These results naturally lead us to a central limit theorem for the Busemann function. Without 
lost of generality, let take a = 57r/4 and let 

where /3 := (cos/?, sin ^). Notice that 

EB{j3,t) = (cos/3 + sin/3)i. 

If /3 G [0,7r/2] then 

B{^,t) = L,m (5.3.4) 
(as process in t). This is just a consequence of (I3.1.ip . On the other hand, if /? G [7r/2,7r] then 

B{^,t) = Y{t)-X{t), (5.3.5) 

where Y{t) and X{t) are two independent Poisson process with intensity 1. 

To obtain ()5.3.5p . note that, by additivity and anti-symmetry (Proposition I1.4T]) . 

B0,t) = Ba{0,{tcosP,0)) + Bo,{{tcosf3,0),tl3) 
= Ba{{t cos 13, 0),t(3) -5c,((tcos/3,0),0) 
= Li((tcos/3,0),t/3) - Li{{t cos l3, 0),0). 

By Theorem [3231 (Burke's Theorem), both t ^ Li(0, (t,0)) and t ^ Li(0, (0,t)) are Poisson 
processes with intensity 1, and t i— )• Li^^{0,t/3) has independent increments. This last statement 
can be seen in the following way: consider the increments Li(0,s/3) and Li{s/3,t/3), for < s < t; 
see Figure 15.11 Then both increments depend on the independent Poisson processes Li restricted 
to (the line segments) 0^4 and AB, and on the Poisson process in the triangle OAB. Burke's 
Theorem for Li now implies that Li restricted to CE is independent of Li restricted to ED. Since 
Li restricted to BE depends on Li on CE, CB and the Poisson process in ECB, whereas Li 
restricted to OE depends on Li on OD, DE and the Poisson process in ODE, and since these six 
processes are all independent, we proved that the increments are independent. 



Corollary 5.3.3 For f3 £ [0, vr] 



Var5(/3,t) 
lim = cos /? — sin m . 

t— s>oo t 
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t cos P 



Figure 5.1: Independent increments for Li. 



Further, for /3 £ [0, 7r/4) 



f B(l3,t) - (cos (5 + sin I3)t \ 
t^oo \ Vcos/^-sin/^ ti/2 - V - 



and for (3 G (7r/4, vr] 



lim 



B{f3,t) - (cos/3 + sin/3)t 



< n = P(iV < n) , 



*^oo \ (-y/sin/3 — cos/3)ti 
where N denotes a standard Gaussian random variable. 

Proof: It follows by combining ()5.3.4p and ()5.3.5p together with Corollary 15.3.21 



66 



CHAPTERS. FLUCTUATION RESULTS 



5.4 Cube-Root Asymptotics for Lx 

In this section we will show how, in the classical Hammersley process, we can use the results we 
have obtained in the previous sections to prove cube-root fluctuations of Lx and L. 

In p] results on the limiting distribution of L were derived, and building on these methods, Baik 
and Rains [3] derived limiting distributions for Lx, using exact formulas for this length, and heavy 
analytic methods to analyze asymptotics of these exact formulas. We were not able to obtain these 
limiting distributions using our methods, but we were able to find the correct scaling. Furthermore, 
we were able to link these results to transversal fluctuations of the longest path. Also, our methods 
have since been extended to other models, where the analytic methods were not (yet) available. 

5.4.1 Upper Bounds 

We start by bounding the fluctuations of Li{t,t). We wish to use Chebyshev's inequality, so the 
following lemma is useful to us: 

Lemma 5.4.1 For all X > 0, we have 

Y&r{Lx{t,t)) < t\X - h + Var(Li(t,t)). 
A 

Proof: Using reflection in the diagonal, we can see that 

Lx{t,t) = Lyx{t,t). 
This means that it is enough to prove the lemma for < A < 1. 

Now couple vx and z^i, so vi contains more points. Suppose for x S ffi and t > that Zx{x, t) > 
Zi{x,t). Then 

L{{Ziix, t), 0), (x, t)) - L{{Zx{x, t), 0), (x, t)) 

> ui{Zx{x,t)) -ui{Zi{x,t)) 

> vx{Zx{x,t)) - vx{Zi{x,t)) . 

This would imply that 

L((Zi(x, t), 0), (x, t)) + ux{Zi{x, t)) > Lx{x, t), 

which is absurd. Therefore, for all x G M and i > 0, we have that Zx{x,t) < Zi{x,t). Using this 
and Theorem 15. 2. H we get 

Yav{Lx{t,t)) = t(i-A) + 2AE(Z+(t,t)) 

< t(i-A) + 2E(Z+(t,t)) 
A 

= t(y-A)+Var(Li(t,t)) 



□ 
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Theorem 5.4.2 There exists a constant ci > such that for all t > and all r > 0, we have 

P(Z,(U)>rEZf(M))<(^^$^(i + i). 

Proof: Note that for any A > 1 and < u < t, 

F{Zi{t,t)>u) = f(3 z > u : ui{z) + L{{z,0),{t,t)) = Li{t,t)^ 
< f(3 z>u:i^i{z) + Lx{t,t)-Lx{z,0) 

>Li(t,t)) 

= ¥(3 z>u:ui{z)-ux{z) 

> Li{t,t) - Lx{t,t)y 

For the inequahty we used the fact that for any A > and any < x < y G M^, we have 

LA(y) >LA(x)+L(x,y), 

and hence, 

LA(y)-^A(x) >L(x,y). 

Note that this is true for any choice of vx^- Since we took A > 1, we can choose I'x as a thickening 
of ui , so that 

Dx-iiz) := i^x{z) - vi{z) 
is in itself a Poisson process with intensity A — 1. This means that 

P(Zi(t,t) >u) <F{i>x-i{u) < Lx(,t,t) - Li{t,t)). 

To have a useful bound for all < n < |t, we choose A such that 



Ei)x-i{u) -E{Lx{t,t) - Li{t,t)) = (A-l)n-t (^A + i 



^ 2 

is maximal. This means that we choose 

Xu = {i-u/tr'/^. 

Some useful elementary inequalities, that hold for all < u < |t, are 

A„ < 2 

EDx^_,{u) -E{Lx^{t,t) - L,(t,t)) > \u^/t (5.4.1) 
A„ - l/\u < 2u/t. 
Using Theorem 15.2. l) Lemma 15.4.11 and (|5.4.ip , we get 

Var(LAjt,t) -Li(t,t)) < 2 {Var(LAjt, t)) + Var(Li(t, t)) } 

< mz+{t,t) + 2t{K-i/\u) 

< 8EZ+(t,t) +4tt. 
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Now we can use Chebyshev's inequality: 



{Zi{t,t)>u) < F{i)x^_i{u)<LxAt,t)-Li{t,t)) 

< P(PA„-iH<EPA„-i(n)-^xV(8t)) 
+ P(LAjt,t)-Li(t,t) >EPA„-i(n)-txV(8t)) 
64t2(A„ - l)u 



< 
+ 



»(LAji,t)-Li(f,t) 



>E{LxAt,t)-L,it,t)) + -) 



64^2 64f2(8EZ+(t, t) + Au) 

— s ' 4 

< 512|4 + ^5SM 



If t > u > |t, we see that 



P(Zi(t,t) > n) < P(Zi(t,t) > -t) < ^ + -1 IA±± 



where we use (j5.4.2|) and choose ci = (4/3)^ x 512. This means that (|5.4.2|) is true for all < u < t, 
if we replace 512 by ci. Since Zi(t,t) < t with probability one, we have that ()5.4.2p is true for all 
ti > 0. The theorem now follows from choosing u = r¥,Z^{t,t). 

□ 



Corollary 5.4.3 

Var(Li(t,t)) EZ+(t,t) 
lim sup iT-jr = 2 hm sup ^—p: < +oo . 

Proof: By Theorem EMI 

poo 

1=1 P(Z+(t,t) > rEZ+(t,t))(ir 

JO 



and hence, 



< 



1/2 

P(Z+(t,t) > rEZ+(t,t))dr 



poo 

+ / P(Z+(t,t) > rEZ+(t,t))dr 

71/2 



'1/2 

1 { f I 1 \ , 1 



< 7^ + ci < / ^ + ^ M'' 



2 1A/2V^' r^J j (EZ+(i,t))3' 
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By Theorem 15.2. 11 this shows Corollary 15.4.31 

□ 

As a consequence, we get: 
Corollary 5.4.4 There exists a constant C2 > such that for all r > 0, 

P(Zi(t,t) >rf2/3) < ^. 

Proof: This follows from Theorem 15.4.21 and the previous corollary. 

□ 



5.4.2 Lower Bounds 

Using the upper bound results from the previous section, we can also get a lower bound result of 
the same order. We will use the sources and sinks representation of the Last Passage Percolation 
model: for x,t > let ux be defined on [—t, x] by 

"^^'^ \ ulit) forzG(0,t]. 

Remember that z^^ is the Poisson process of "sinks" on the t-axis, with intensity 1/A. Thus I'x, 
restricted to [0,x], is a Poisson process with intensity A (sources) while I'x restricted to (0, i] is a 
Poisson process with intensity 1/A (sinks). 

Denote z := (z,0) for z G [0,x] and z := (0, z) for z G [-t,0). By (!2.3.2p we have that 

Lx{x,t)= sup {vx{z) + L{'L,{x,t))} . 

z&[-t,x\ 

Define 

Zx{x, t) := sup {z G [-t, x] : i?x{z) + L(z, (x, t) = Lx{x, t)} , (5.4.2) 

and 

Z'^{x, t) := inf {z G [-t, x] : ux{z) + L(z, (x, t) = La(x, t)} . (5.4.3) 
Thus, Zx and Z'-^ are exit points for i'x- 

By definition, and using the invariance under reflection, we can see that 

Z'^ix,t) < Zxix,t) and Zx{x,t) = -Z[^^{t, x) . (5.4.4) 

It also clear that 

Z+{x,t) = Z+{x,t). 

By using the same reasoning as in the proof of Proposition 12.3. 1| the following local comparison 
property holds. For any X2 > xi > x, if Zx{x,t) > then 

L(0, ix2,t)) - L(0, (xi, t)) < Lx{x2,t) - Lxixi,t) . 
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Furthermore, for any xi < X2 < 2; , if Z'^{x,t) < then 

L{0,ix2,t)) - L{0,{xi,t)) > Lxix2,t) - Lxixut). 

From now on we denote 

A{z,t) ■.= L{0,{t,t))-L{0,{t-z,t)), 

so that, if Zx{t,t) > then 

A{z,t)) <Lx{t,t)-Lx{t-z,t), (5.4.5) 

and if Z'^{x,t) < then 

A{z,t)>Lx(,t,t)-Lx{t-z,t). (5.4.6) 

Lemma 5.4.5 For r > let 

r r 
A+ = 1 + — 7^- > 1 and A_ = 1 -7^ < 1 . 

Then there exist constants ci , C2 > such that 

F{Zx4t,t)<0) < ^, 

and 

p(zL(t,t)>o) <^, 

for all sufficiently large t and r G [l,t^/^). 

Proof: By rescahng and translation invariance, we see that 

Zx{t,t) = ^Zi{Xt,t/X) and Zx{x + h,t) = Zx{x,t) + h . 
A 

So, by Corollary 15.4.41 

F{Zx4t,t)<0) = P(Zi(A+t,i/A+) <0) 
< F{Z[{X+t,t/X+) <0) 
= F(Zi(t/A+,A+t) > 0) 
= P(Zi(A+t,A+t) > (A+-l/A+)t) 

'Zi{X+t,X+t) > 2^2/3 - 
= 0(r-3), 

for large enough t and r G [l,t^^^). Analogously, 

¥ (^Z'x_{t,t) > Oj = F{Z[{X.t,t/X-) > 0) 

< P(Zi(A_t, t/A-) > 0) 

= F(Zi(t/A_,t/A_) > t(l/A_ - A_)) 

= 0(r-3), 

for large enough t and r G [l,t^^^). 

□ 
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Lemma 5.4.6 Fix 6 > 0. Then 



limsupP sup {ui{z) - A(z, t)} > 5t^/^ = 0{e^) . (5.4.7) 

*^°<3 \2e[0,et2/3] J 

Proof: Let A = 1 — rt~^^^, for some r > 0. We introduce a Poisson process i^x of sources on the 
X-axis, independent of vi, with intensity A and denote 

iyi{z) = Lx{t,t)-Lx{t-z,t). 

Clearly, vj^ is also a Poisson process of intensity A, independent of i^i. Using (I5.4.6|) . we get that 
the probability in the left hand side of ()5.4.7p is bounded by 



sup {ui{z) - uiiz)} > St'/' +F{Z'^{t,t)>0) , 

,0<2<et2/3 / 



which is equal to 



P sup {i^i{z) - ux{z)} > 6t'/' + P {Z'^{t, t)>0) . (5.4.8) 

We will treat these two terms separately. We start with the second term. By Lemma 15.4.51 

P(Z^(t,t) >0) =0(r-3) 

for large enough t and r € [l,t'/'). 

The first term in (j5.4.8p concerns the hitting time for the difference of two independent Poisson 
processes. The process 

{z^r'/'[Mz)-i.x{z)]} 

converges, as t — )• oo (in the topology of uniform convergence on compacta), to the drifting Brownian 
motion process, namely 

Wr{z) := W{2z) + rz, 
where W is the standard Brownian motion on M+. We now get, for r < 5/e, 

P ( sup Wriz) > 5] < P ( sup W{z) > 5-re] 

\ze[0,e] J \ze[0,2e] J 

^( W ^ 5-re\ 

= P sup 77^\/2^ > 

I 6 — re \ 

= P sup W{z) > 



e[o,i] V2e 
e-^"/^du. 



2 



V2s 
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CO poo 

e-'^ /^du = 2 / -^du 



Taking r = 5/{2e), we get 



TT J s-rt J s \/2ir 

2^/8ie-'^'/(i6^) 

using Mills' ratio approximation for the tail of a normal distribution in the last step. Using Donsker's 
Theorem, this means that, with this choice of r, our estimate for the second term in ()5.4.8|) is 
dominant. 

□ 



Lemma 5.4.7 Let r] > 0. Then there exist 6 = 6{t]) > and e = e{6) > such that 

limsupP sup {i^i{z) - A(z, t)} < 5t^/^ < r? . (5.4.9) 

i^oo y2>et2/3 J 

Proof: The proof of this lemma is similar to the previous one, but now we have a different choice 
of A and r. As before, we introduce a Poisson process u\ of sources on the x-axis, independent of 
i/i, with intensity A = 1 + rt~^/^ and denote 

u{{z) = Lx{t,t) - Lx{t - z,t) . 

Clearly, z^^ is also a Poisson process of intensity A, independent of ui. Using (I5.4.5p . we can 
bound the left hand side of (j5.4.9p by 



sup {vi{z) - vi{z)] < 6t^/^ +F{Zx{t,t) < 0) 

ze(ei2/3,ri2/3] / 



that is equal to 



P sup {iJi{z) -ux{z)} <6t^/^] +F{Zx{t,t) <0) . (5.4.10) 

\2e{ei2/3,rt2/3] / 

(Notice that we also have restricted z to (et'^/^ ,rt'^/^]). 
By Lemma 15.4.51 we see that 

P(ZA(t,t) <0) =0(r-3), 

for large enough t and r G [1, t^^^). This means that we can take r big enough, such that the second 
term in ()5.4.10p is smaller than r//2. 
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Now consider the first term in (|5.4.10p . After rescaling, the term can be written as 

P ( sup {i^i{z) - ux{z)} < I 

\2e(et2/3,rt2/3] / 

= P ( sup {t-^/^[ui{zt^/^) - i^xizt^^^)]] < 5] . 

We know that the process 

{z^t-'/^[u,{zt'/')-ui{zf/')]} 

converges, as t — >■ c«, to a drifting Brownian motion process, namely 

Wriz) := W{2z) -rz ,z>0, 

where W is standard Brownian motion on M+. A standard application of Donsker's Theorem then 
shows that 



lim P sup {ui{z) - u{{z)} < 6t^/^ 

= fI sup Wriz) < 5] 
\z£{e,r] J 

< P ( sup Wr{z) <6] + fI sup Wr{z) > 5] . 
\z<a[0,r] J \^G[0,£] J 

Remember that we have fixed r > already. Now we choose 6 = 6{r], r) small enough, such 
that 

P ( sup Wr{z) <6] < r?/4. 

\ze[0,r] J 

For fixed r and 6{rj,r), we can now choose e > small enough (as in the previous estimates), such 
that 

P ( sup Wriz) > 6] <r]/A. 



□ 



Now we are able to prove the main result of this section: 
Theorem 5.4.8 



lim lim sup P (O < < et^^^ 

Proof: Choose rj > 0. It is enough to find e > such that 

limsupP (o < Ziit,t) < et^/^) < r] . (5.4.11) 
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The probability in the left hand side of (|5.4.1ip is bounded by 



sup {ui{z)+L{{z,0),{t,t)))} 



< sup {z^i(z) + L((z,0),(t,t)) 

0<^<ei2/3 



For any given (5 > 0, this is bounded by 



sup {vi{z) - A(z,t)} < 5t^/^ 



+ P( sup {z^i(z) - A(z,t)} > . 

\0<^<ei2/3 J 



We note that the random function •) on M x (0,oo) is independent of vi. Furthermore, 



V 



{z ^ L(0, {t, t)) - L{{z, 0), {t, t))} = {z^ L{0, {t, t)) - L(0, {t - z, t))} 



(as processes). Therefore, (|5.4.1ip follows by combining Lemma 15.4.71 together with Lemma 15.4.61 

□ 



Corollary 5.4.9 



liminf'^^(L"'=liminf^"'^,'L"»>0. 



Proof: Since, with probability 1, either Zi{t,t) > or Zi{t,t) < 0, we have that 

1 = > 0) + P(Zi(t,t) < 0) 

< P(Zi (t, t)>0)+F{Z[ {t, t) < 0) 
= 2P(Zi(t,t) > 0) , 



and hence, 



which shows that 



< P(Zi(t,t)>0) 
= P(Zi(t,t)>0) 

< Zi (t , t) < et^/s] + p (zi {t,t)> ei2/3 



< < £^2/3 



< 



EZ+(t,t) 



t2/3 



By Theorem 15.4.81 one can choose Eq so that 

p(o<Zi(t,i) <eot'/') < ^, 
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for large enough t, which imphes that 



eo ^ EZt{t,t) 



4 t2/3 ' 

for large enough t. Theorem 15.2.11 completes the proof of the corollary. 



□ 



5.5 Cube- Root Asymptotics for L 

In the preceding sections it was shown that, for the classical Hammersley model in equilibrium, the 
variance of of Li{t,t) is of order t^/^. The key result to extend the cube-root asymptotics to L is 
the following theorem: 

Theorem 5.5.1 There exists a constant ci > such that for all r > 

p(|Li(t,t)-L(t,i)|>rti/3)<-|^. 

Before proving Theorem 15 . 5 . 1 1 we have: 
Corollary 5.5.2 There exists a constant C2 > such that 

E\L{t,t) -2t\ < C2t^/^ . 

Proof: Indeed, 

¥.\L{t,t) -2t\ < E|L(t,t) - Li(t,t)| +E|Li(t,t) - 2t| 

< + r-s/^drj + ^VarLi(t,t) . 

By Corollary 15.4.31 this proves Corollary 15.5. 2[ 

□ 

To prove Theorem 15.5.11 we require the following lemma, the proof of which uses again the local 
comparison argument. 

Lemma 5.5.3 There exist constants ci, C2 > such that for all t > and all r G [1, cit^/^] 

P( sup {u^{z) - A{z,t)} > rt'/^) < ^ . (5.5.1) 

Proof: First, let K = r^/^^. As in the proof of Lemma l5.4.6t for A = 1 — ut^^^^ we get that the 
left hand side of (|5.5.ip is bounded by 

pf sup {ui{z)-ux{z)} >rt'/A +F{Z'^{t,t)>0) , 

\\z\<Kt2/3 J 
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and 

F{Z[{t,t)>0) =0{u~^), 

uniformly for u G [l,t^/^/2]. To deal with the first term, we first note that w.l.g. we can restrict z 
to [0,i^t2/3] (foj. ^ g the argument is similar). Also that 

z G [0, K] ^ Mu{z) := r^^^ {uiizt^^^) - uxizt^^^)} - uz , 
is a zero mean martingale, and that 

E{Mu{zf) < z + (1 - X)z < 2z . 

Hence, 

P( sup {ui{z) -uxiz)} > rt^/^] 

= P I sup {Mu{z) +uz} >r] 
\zelO,K] J 

<P( sup {Mu{z)} >r-uK] . 
\zelo,K] J 

Taking u = r/{2K) and using Doob's sub-martingale inequality, we get 
P( sup {iJi{z) -ux{z)} > rt^/^] 

< P f sup {Mu{z)} > r/2 ) 
\zelo,K] J 

= P( sup {M„(z)2} >rV4) 
\zelo,K] J 



^2 ^2 ^4/3 

This means that, for u ~ rK~^ = r^/^"^ g [l,t^/'^/2], our estimate for the second term is dominant, 
and the lemma is proved. 

□ 

Proof of Theorem [5311} By Corollary [533] 

p(Li(t,t)/ sup {i/i(z) + L((0,z),(t,t))} I 

\ |2|<Xt2/3 J 

< 2P (Zi{t,t) > Kf/A = 0{K-^) . 
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So taking K = r^/^^ and using Lemma 15.5.31 (and translation invariance) we have, for all r E 

p(|Li(t,t)-L(0,(t,t))| >rti/3^ 

^Li{t,t) - L{0,{t,t)) > rt^/^ 



sup {ui{z) - A{z,t)} > rt^l^ 

|z|<r7/12j2/3 ^ 



+ 2¥(Zi{t,t) > Kt^/^ 

we note that 

P(^Li(t,i)-L(0,(t,t)) >rt^/3) < P(^Li(t,t) > 

/ 5/6 \ 

where Pt is a Poisson random variable with expectation t. This implies that 

'Li{t,t) -L{0,{t,t)) > rt^/'^^ 



tends to zero faster than any negative power of r, if r > t^/^ ^ uniformly in all large t. Hence, we 
can conclude the proof of Theorem 15.5.11 

□ 
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